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ABSTRACT 

V 

This  report  summarizes  the  results  of  a  series  of  studies  to 
establish  engineering  methods  for  estimating  the  flow  fields  and  aero-,,  ,, 
dynamic  forces  over  sharp  cones  during  hypersonic  re-entries^'  The  effects 
of  ablation  and  turbulent  flow  are  not  considered. 

The  flow  fields  and  aerodynamic  forces  on  sharp  cones  at  zero 
angle  of  attack  are  first  examined  by  a  comparison  of  various  techniques 
of  determining  the  boundary  layer  displacement  thickness  and  shock  wave 
geometry.  Good  agreement  is  obtained  between  a  modified  viscous  interaction 
theory  and  experimental  force  data.  The  work  is  then  extended  using  vis¬ 
cous  interaction  theory  to  obtain  a  correlation  of  experimental  data  for 
flow  conditions  ranging  from  cdntinuum  to  free  molecule  flow.  The  approxi¬ 
mate  form  of  the  interaction  theory  derived  shows  that  the  parameter, 
is  the  proper  correlating  parameter  for  the  drag  of  slender 
cones  over  a  wide  range  of  Mach  numbers  and  Reynolds  numbers.  Excellent 
agreement  is  obtained  for  experimental  data  correlated  according  to  this 
parameter . 

/  The  zero  angle  of  attack  case  is  then  extended  to  include  the 
effects 'of  the  angle  of  attack  range,  0  to  90  degrees,  on  the  aerodynamic 
forces  and  the  flow  properties  at  the  most  windward  streamline.  Comparison 
was  again  made  between  the  theory  and  experiment,  where  generally  good 
agreement  was  found  to  exist  . 

’ft 

This  work  was  conduce^  by  Aeronutronic  Division  of  Ford  Motor 
Company  under  Contract  AF  04(6941-23  to  the  Ballistic  Systems  Division 
of  the  United  States  Air  Force. 
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CHAPTER  1 

INTRODUCTION 


The  development  of  improved  heat  shields  and  a  desire  for  in¬ 
creased  terminal  velocities  have  led  to  increased  interest  in  slender 
body  configurations  for  re-entry  vehicles.  The  investigation  of  the  flow 
fields  and  aerodynamic  characteristics  in  the  high  speed,  high  altitude 
regime  is  complicated  by  the  effects  of  viscous  interaction.  The  con¬ 
siderable  distortion  of  the  flow  field  over  slender  bodies  in  hypersonic, 
low-density  flow  due  to  boundary  layer  displacement  effects  has  been  ex¬ 
tensively  examined  qualitatively1.  Quantitative  results  have  been  presented 
for  special  cases,  but  numerical  answers  directly  applicable  to  the  cone 
under  re-entry  conditions  are  scarce.  Since  it  is  not  generally  possible 
to  match  wind  tunnel  Mach  number  and  Reynolds  number  test  conditions  to 
those  encountered  in  flight,  there  has  been  a  heavy  emphasis  on  the  use 
of  theory  for  predicting  flight  characteristics.  In  addition,  hypersonic 
theory  has  not  yet  been  developed  to  the  point  where  solutions  are  avail¬ 
able  for  all  density  regimes  from  free  molecule  to  continuum  flow. 


In  order  that  engineering  estimates  might  be  made  of  flow  field 
and  aerodynamic  characteristics,  without  resorting  to  high  speed  computing 
machines,  it  is  common  to  make  certain  assumptions  to  simplify  the  complex 
equations  of  motion.  For  example,  the  use  of  constant  density  theory  in 
hypersonic  flow  has  provided  results  for  inviscid  conical  flow^  which  is 
highly  accurate,  yet  gives  analytical  expressions  adaptable  to  hand  cal¬ 
culations.  The  assumption  of  a  constant  density  should  in  no  way  be  taken 
to  imply  an  incompressible  flow,  but  rather  to  be  particularly  applicable 
in  hypersonic  flows  where  thin  shock  layers  occur.  The  general  premise 
underlying  the  present  work  is  that  the  local  flow  field  can  be  sufficiently 
described  for  engineering  purposes  if  the  boundary  layer  growth  and  re¬ 
sultant  shock  wave  can  be  computed  independently  of  any  considerations 
other  than  free-stream  conditions  and  body  geometry.  A  considerable  amount 
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of  Information  about  the  flow  field  can  be  specified  from  the  local  slope 
of  the  shock  wave.  This  "inverse"  method  has  two  major  advantages; 

(1)  boundary  conditions  are  provided  with  the  deletion  of  one  independent 
parameter  so  as  to  do  away  with  the  usually  required  simultaneous  solution 
of  several  equations,  and  (2)  the  matching  of  the  inviscid  and  viscous 
flow  regions  is  automatically  satisfied  by  the  use  of  the  equivalent  body 
in  obtaining  the  shock  wave  position.  The  major  disadvantage  is,  of  course, 
the  uncertainty  in  the  basic  assumption  Itself.  In  the  present  case, 
existing  methods  (compressible  laminar  boundary  layer  and/or  viscous  inter¬ 
action  theories)  or  developments  thereof,  will  be  used  to  estimate  the 
boundary  layer  growth. 

This  continuum  flow  solution,  including  the  effects  of  viscous 
interaction,  provides  a  firm  basis  for  the  understanding  of  the  behavior 
and  trends  exhibited  as  the  density  is  decreased.  In  the  absence  of  a 
solution  for  the  complete  flow  regimes,  a  correlation  of  experimental  data 
appears  most  promising  for  engineering  designs.  This  is  also  accomplished 
in  the  present  investigation  for  the  zero  angle  of  attack  case. 
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CHAPTER  2 

REVIEW  OF  VISCOUS  INTERACTION  THEORY 


2 . 1  GENERAL 


At  low  Reynolds  numbers  and/or  high  Mach  numbers  the  boundary 
layer  becomes  significantly  large  that  it  produces  an  outward  deflection 
of  the  stream  surfaces  and  perturbs  the  local  inviscid  flow  field.  This 
perturbation  of  the  Inviscid  flow  field  results  in  an  "induced"  pressure 
change  on  the  body  which,  in  turn,  affects  the  growth  of  the  viscous  boundary 
layer  and  also  the  local  skin  friction.  This  problem  of  boundary- layer- 
shock  wave  interactions  has  been  under  intensive  investigation  for  some 
time  and  a  number  of  solutions  of  varying  degrees  of  accuracy  are  available. 
Because  of  the  mutual  effects  of  the  boundary  layer  on  the  shock  wave  and 
the  shock  wave  on  the  boundary  layer,  it  is  to  be  expected  that  the  most 
accurate  interaction  solutions  are  those  of  the  type  involving  an  iterative 
procedure.  These  procedures  are  cumbersome  and  time  consuming,  therefore 
a  number  of  approximate  theories  have  been  developed  which  are  applicable 
over  certain  ranges.  The  original  theory  developed  by  Lees  and  ProbsteinJ 
is  used  as  the  basis  for  these  interaction  theories.  It  has  been  amply 
examined  (Reference  4-6)  ,  and  only  a  brief  review  of  the  applicable  equa¬ 
tions  will  be  given  nere.  An  iterative  procedure  will  be  discussed 
subsequently. 


2.2  SHARP  FLAT  PLATE 

Before  proceeding  to  the  cone,  the  case  of  the  sharp  flat  plate 
will  be  examined.  In  the  case  of  the  flat  plate,  it  has  been  found  con¬ 
venient  to  divide  the  interaction  effects  into  two  asymptotic  regions 
called  the  weak  and  strong  interaction  zones. 
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For  the  weak  interaction  zone  the  effects  produced  by  the  self 
induced  pressure  gradient  are  considered  to  be  essentially  perturbations 
superimposed  on  an  already  existing  uniform  flow.  Under  these  conditions, 
the  boundary  layer  growth  is  essentially  unaffectedly  the  pressure  gradient. 
For  the  strong  interaction  region,  on  the  other  hand,  the  increment  in  flow 
deflection  caused  by  the  boundary  layer  is  of  the  same  order  of  magnitude  as 
that  caused  by  the  body  and  the  pressure  gradient  is  considered  to  have  a 
direct  influence. 

The  usual  division  between  the  strong  and  weak  Interaction  regimes 
is  written: 

a)  x  <<1  weak  interaction 

b)  x  >  >!  strong  interaction 

where  x^s  t^e  hypersonic  interaction  parameter  defined  as: 

X-  ^r- 

v  Re 

where  C  is  the  coefficient  in  the  linear  viscosity-temperature  relationship. 
This  relationship  is  used  to  simplify  the  boundary  layer  equations  to  obtain 
closed  form  solutions  for  the  boundary  layer  parameters.  Because  temperature 
affects  both  the  density  and  viscosity  of  a  fluid,  it  is  apparent  that  the 
drag  on  a  body  (from  interaction  and  skin  friction  sources)  is  a  function  of 
the  temperature  distribution  through  the  boundary  layer.  This  temperature 
distribution  is  related  to  the  body  wall  temperature.  While  exact  numerical 
solutions  for  the  boundary  layer  characteristics  for  various  temperature 
distributions  are  available^,  it  has  been  found  convenient  from  an  analytic 
point  of  view  to  use  a  linear  variation  of  viscosity  with  temperature  of 
the  form: 


^ref  Tref 

M  T 

It  was  originally  proposed*  that  the  reference  conditions  be 
evaluated  at  the  body  wall.  Now,  since  the  temperature  (and  viscosity)  vary 
by  large  amounts  across  the  boundary  layer,  it  is  apparent  that  a  considerable 
variation  in  C  may  also  occur.  It  would  therefore  seem  that  the  most  realistic 
evaluation  of  C  would  be  at  some  mean  boundary  layer  temperature,  rather  than 
the  wall  temperature. 


2-2 


AERONUTRONIC  OIVISION 


A  number  of  authors  have  submitted  various  recommendations  as  to 
how  the  reference  conditions  should  be  evaluated. 

g 

Talbot,  Koga  and  Sherman  use: 

T  ,  -  T'  =  T  -  0.486  Pr1/3  (T  -  T.  )  -  0.273  Pr  (-^ r )  M2  T 
ref  w  w  Aw'  2  c  c 

where  Pr  ■  Prandtl  number 

TAw  *  Adiabatic  wall  temperature 

Cheng,  Hall,  Golian,  and  Her tzberg?  recommend  simply  the  average 
temperature  across  the  hypersonic  boundary  layer: 


Trcf  ■  (1  +  3  (Tu/T0)  J  / 6 

T0  '  Fr”  *"*"  temperature 

Sommer  and  Short ^  give: 

T  ,  ■  T"  -  T  +  0.032  T  M2  +  0.58  (T  -  T  ) 
ref  c  c  c  w  c 

Eckert’s  semi  empirical  reference  enthalpy*^  gives: 

T  .»!*:■  0.5  (T  +  T  )  +  0. 22  Pr1 /2  (-^)  M2  T 
ref  '  w  c  2  c  c 

It  appears,  on  the  basis  of  experimental  evidence,  that  Eckert's 
reference  enthalpy  method  gives  the  be^t  results.  However,  the  numerical 
differences  between  the  various  reference  temperatures  are  generally  quite 
small.  Table  I  compares  reference  temperature  calculated  by  the  above 
four  methods  for  various  assumed  Mach  numbers  and  Prandtl  numbers. 

For  the  remainder  of  this  discussion  only  Eckert's  reference 
temperature  method  will  be  considered.  In  the  case  of  a  dissociated 
boundary  layer  flow,  a  reference  enthalpy  will  be  used, 
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h*  -  0.5 


<h  + 

w 


V 


1/2 

+  0.22  Pr  '  U 


72 


where  the  reference  temperature,  T* ,  Is  to  be  obtained  from  a  Mollier 
diagram  at  (h*,  P2) . 

Once  the  reference  temperature  is  established,  the  viscosity 
proportionality  constant  Cc  may  be  evaluated  by  means  of  either  the  power 
law  or  the  Southerland  formulation.  The  power  law  leads  to  a  simpler 
algebraic  relation,  but  its  applicability  is  limited  to  more  restricteo 
temperature  range  than  the  Southerland  equation. 

From  the  power  law  (u>*  0.76): 

(applicable  in  the  range  170°K  < T  <500°K) 


whereas  for  the  Southerland  formula: 
(100°K  <  T  <1900°K) 


or 


C 

c 


C 


c 


I  + 


110°K  Tref 


1  + 


ref 

110°K 


ref 


ref 


ref 


1/2 


T  +  199 
c _ 

T  ,  +  199 
ref 


1/2 


,  T  a  °K 

T  a  °r 


It  can  be  shown  that  the  use  of  Eckert's  equation  to  determine 
the  reference  temperature  together  with  the  Southerland  formulation  to 
determine  the  viscosity  coefficient  gives  boundary  layer  solutions  that 
are  within  about  4  percent  of  the  correct  values  calculated  numerically 
by  Van  Driest. 
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The  Induced  pressures  produced  on  a  flat  plate  or  wedge  by  the 
phenomenon  of  viscous  interaction  has  received  extensive  examination.  For 
completeness,  the  results  presented  by  Kayes  and  Probstein  will  be  re¬ 
produced  here.  The  subscript,  0,  refers  to  the  original  inviscid  value. 
The  induced  pressure  on  a  wedge  in  the  weak  interaction  region,  including 
heat  transfer,  is  (to  second  order): 


4-  -  1  ♦  7J,  d„  ♦  UX±U 


1  0  M0  A0 


V 


where,  _ 


XQ  -  Mj  a/c-  /  */Re 


-  -4-  ~  +  (  7-1)  B=  A  7L+(7-l)B 


0 


”0 


{ 


1  +7  K 


V<T)2  +  +  ^]) 


-1 


0±l  K  +  1  +  2  Kb 
2  b  Vl+  <W  K?  J 


2.1 


2.2 


2.3 


2.4 


m  y, 


2.5 


For  the  flat  plate  at  zero  angle  of  attack,  K,  -  0,  then  J.  -  J„  *  1  and 
Equation  2.1  reduces  to:  1  Z 


P 

P 


0 


1  +  G 


+ 


7+1 
4  7 


G 


2.6 
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where,  T 

G  -  0.166  7(7  -1) (2.605  ^  +  1),  Pr  -  1  2.7 


The  constants  A  and  B  in  Equation  2.2  are  given  in  the  following  table: 


Pr 

A 

B 

1.0 

0.725 

0.865 

0.968 

0.166 

0.145 

For  Pr  -  0.725, 

T 

G  -  0.145  7(7  -1)(3. 24  rr  +  1)  2.7a 

Tt 

For  a  wedge  withKb»l  (say  Mw0b  >5),  Jp  S  2/7  Kb,  J2  ~ 

4/ 7  (7  +  1)  Kb ,  and  0b=-^5/  y  (  V-l)  .  The  functions  Jp  and  J2  are 
plotted  versus  Kb  in  Figure  2-1. 


In  the  region  of  strong  interaction  (to  first  order)  and  at 
zero  angle  of  attack, 


P 

P 


0 


Xn  +  0.76 


2.3 


The  constant,  0.76,  was  derived  for  the  case  of  an  insulated  plate,  and  is 

assumed  to  hold  approximately  for  all  values  of  Tw/Tt. 

T 

w 

The  function(— )  is  given  in  th°  following  table  for  Pr  =  1 : 


T 

w 

Tt 

J 

0 

y  *  1  .4 

1.67 

0 

0.149 

0.261 

0.2 

0.232 

0.403 

0.6 

0.377 

0.696 

1.0 

0.514 

0.921 
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For  Pr  *  0.7,  calculations  have  been  made  which  give 

$ p<^>  -  -§  &  4  7(7+  1)  (7  -  1)  \ 


The  integral,  ,  is  presented  in  Figure  2-2  as  a  function  of  Tw/Tt. 

The  induced  pressures  for  the  two  regions  of  interaction  are 
presented  in  Figure  2-3  to  various  orders.  The  closest  matching  solutions 
is  the  second  order  weak  interaction  and  first  order  strong  interaction. 
This  is  given  in  Figure  2-4.  It  is  noted,  however,  that  at  the  largest 
values  of  ^  ,  the  effects  of  slip  are  ignored,  which  would  reduce  the 
induced  pressures. 


The  viscous  interaction  theory  for  the  boundary  layer  dis¬ 
placement  thickness  is  presented  in  Reference  1  as 
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Pr  ■=  0.725 


where 


X  «1 


T*  ■  <VHt>  X  /Morig,  Pr  -  liX  »1 
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where  1 6 (Hy/H^)  is  a  function  of  the  wall  temperature  and  is  given  in  the 
following  table  for  the  zero  order  theory: 


Hw 

Ht 

<VHt> 

7  -  1.4 

7  -  1.67 

0 

0.397 

0.457 

0.495 

0.568 

mEM ■ 

0.632 

0.726 

BSm 

0.738 

0.858 
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FIGURE  2—2 .  INTEGRAL  I  FOR  SIMILAR  SOLUTION 
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FIGURE  2-3.  COMPARISON  OF  INTERACTION  SOLUTIONS  FOR  INDUCED  PRESSURES  ON  A 

SHARP  FIAT  PLATE 
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For  Pr  -  0.7,  Reference  6  replaces  Equation  2.11  for  the  strong 
interaction  region 

^  -  #|  v,)  X  1/2 /",ne  2 

where 


^(W 


4  o- 1)  j5 

J  n/7  (  7+  1) 


I 


1 


Equation  2.10  is  not  an  interaction  solution  in  the  sense  that  the  weak 
induced  pressures  are  accounted  for;  but  it  is  derived  from  the  zero 
pressure  gradient  solutions  of  Crocco,  with  a  linear  viscosity-temperature 
law  and  constant  Prandt  number.  Hayes  and  Probstein  have  shown  that  the 
induced  pressure  has  no  first-order  effect  on  the  slope  of  the  boundary 
layer,  but  produces  only  a  uniform  reduction  in  thickness.  This  reduction 
can  be  calculated  from  results  presented  by  Lees  and  Probstein,  providing 
the  first-order  solutions  (in  the  assumed  velocity  and  temperature  series 
expansions)  to  the  energy  and  momentum  equations  are  known.  These  solutions, 
as  obtained  by  Lees  and  Probstein,  have  been  utilized  to  estimate  a  correction 
to  the  boundary  layer  thickness,  6*(°) ,  as  expressed  by  Equation  2.10.  To 
the  accuracy  of  the  hand  calculations  (including  several  graphical  integrations), 
the  corrections  for  Pr  =  1  may  be  expressed  as  follows: 

Tw/Te  *=  0  (highly  cooled  body) : 

6*(1)*  6*(°>  (i  -  0.01  X  Q) 

Tw/Te  =  TAw/Te  (insulated  wall):  2.13 


5*(1)  =  6*(0)(1  -  0.10  XQ) 


For  the  weak  interaction  assumptions  (")(q<1),  it  is  seen  that  the  correction 
amounts  to  less  than  one  percent  at  T^/Tg  ■  0,  and  less  than  ten  percent  for 
the  insulated  plate  case. 

The  displacement  thickness  is  presented,  rather  than  the  boundary 
layer  thickness;  but,  in  hypersonic  flow  where  the  densities  in  the  boundary 
layer  are  very  low,  6  and  6*  may  be  taken  to  be  approximately  equal.  This 
is  a  considerable  increase  in  6  *  over  the  incompressible  flow  case  where 
6*  =  6/3.  An  estimate  of  the  deviation  of  6*/6  from  a  value  of  one  may 
be  made  as  a  function  of  Me  and  Tw/Te.  For  a  perfect  gas  and  constant  pres¬ 
sure  across  the  boundary  layer,  the  results  of  the  Crocco  integral  give  the 
density  distribution  in  the  boundary  layer  as: 


2-12 


I 

I 

I 

1 

I 

1 

I 


o 


■tf 

00 

00 


°°d/d 


2-13 


FIGURE  2—4 .  INDUCED  PRESSURES  ON  A  SHARP  FLAT  PLATE  DUE  TO  VISCOUS  INTERACTION 
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From  Che  definition  of  displacement  thickness, 


Then,  for  T  /T  »  1  (local  hypersonic  flow), 
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The  variation  of  8*/8  with  local  Mach  number  for  a  linear  velocity  profile 
and  various  values  of  Tw/Tfi  is  given  in  Figure  2-5.  This  implies  that  in 
the  very  strong  interaction  regime  where  the  local  Mach  number  may  drop 
considerably,  due  to  the  highly  curved  and  very  thick  boundary  layer,  the 
assumption  of  6*/8  a  1  could  be  in  appreciable  error.  The  hypersonic  flow 
criteria  must  be  considered  to  be  local  hypersonic  flow  for  8*  ~  5.  The 
two  thickness  will  be  used  interchangeably  in  the  remainder  of  this  report. 


Various  procedures  have  been  devised  to  provide  a  smooth  transi¬ 
tion  from  the  weak  to  strong  interaction  region.  Approximate  solutions, 
such  as  integral  and  inverse  methods,  can  be  constructed  which  inherently 
include  a  continuous  transition.  The  solutions  given  by  Equations  2.10  to 
2.12  consider  the  limiting  regimes  of  weak  and  strong  interactions.  A 
simple  transition  can  be  achieved  by  matching  the  individual  solutions 
at  their  common  intersection  or  point  of  closest  approach.  This  has  been 
done  by  Bertram. 5  His  results  are  as  follows: 
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FIGURE  2-5.  RATIO  OF  DISPLACEMENT  TO  BOUNDARY  LAYER  THICKNESS  IN  HYPERSONIC  FLOW 
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2.3  SHARP  CONE 

For  the  case  of  a  flat  plate,  a  number  of  noniterative  Interaction 
solutions  exist  both  for  the  weak  and  strong  Interaction  regimes.  For  the 
case  of  a  cone,  however,  only  the  weak  interaction  solution  Is  presently 
available.  It  might  be  supposed  that  this  would  lead  to  some  difficulty  at 
the  lower  Reynolds  numbers,  since  under  these  conditions  hypersonic  theory 
would  predict  a  flow  field  of  the  strong  interaction  type.  However,  it 
will  be  shown  later  that  at  these  lower  Reynolds  numbers,  other  effects 
such  as  slip  flow  already  predominate  to  a  degree  that  interaction  theory 
alone  is  no  longer  valid. 

The  conical  flow  weak  interaction  theory  to  be  reviewed  here  is 
that  due  to  Lees  and  Probstein  and  subsequently  amplified  in  Reference  4. 
Probstein's  approach  is  based  on  an  asymptotic  series  expansion  of  the 
boundary  layer  flow  quantities  in  terms  of  ascending  powers  of  the  hyper¬ 
sonic  interaction  parameter  X-  The  zero  order  equations  and  boundary 
conditions  for  this  series,  (after  applying  the  Mangier  transformation), 
are  identical  to  the  flat  plate  case,  and  each  higher  order  equation  is 
expressible  in  terms  of  coefficients  which  are  functions  only  of  lower- 
order  solutions  in  the  expansion. 

From  these  equations,  a  first-order  solution  to  the  interaction 
effects  on  the  cone  drag  are  obtained  using  the  assumptions  of  Prandtl 
number  unity  and  a  perfect  gas  of  constant  specific  heat. 

According  to  Probstein's  theory,  the  interaction  pressure  is 
assumed  to  be  a  unique  function  of  the  local  total  flow  deflection  (0) 
and  may  be  expressed  as  a  Taylor  series  expansion  in  powers  of  d8*/dx. 

For  both  the  wedge  and  the  cone,  the  expansion  is  in  the  follow¬ 
ing  form: 


p  -  1  +  1 

M(p/pm)| 

id5*\  i  |d2  P/PJ 

|  d5*\2 

2.17 

p  p  /p 

C  C  00 

1  de  /e=eb 

U  1  2  Pc/pJ  de2  *e=e 

D 

Url  +  • 

where  is  the  equivalent  body  angle  and  the  subscript  c  refers  to  the 
original  conditions  associated  with  uniform  inviscid  conical  flow  (Taylor  - 
Maccoll  values) . 

While  the  coefficients  of  d5*/dx  can  always  be  evaluated  numer¬ 
ically  from  the  Taylor  -  Maccoll  theory  it  is  inconvenient  from  the  point 
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12 

of  view  of  obtaining  an  analytic  interaction  theory.  Lees,  using  hyper¬ 
sonic  similarity  theory,  has  obtained  an  algebraic  relation  for  the 
inviscid  surface  pressure  on  an  unyawed  cone  that  is  quite  accurate  within 
the  limitations  that  K  =  0C  >  1  and  6C  <  24°.  This  relation  is 


r?T  «2.  -  *>  +  '  »2 


+  1 


Y  -  1  +  2/k; 


s  J 


+  1 
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where  Ks  =  M*,  0g,  08  is  the  shock  angle,  and  6C  is  the  cone  angle.  The 
relationship  between  Ks  and  K  is  given  by: 


K 


8 


_  1± 1 
V  +  3 


2.19 


The  results  given  for  the  wedge  were  based  on  a  solution  for  the  coefficients, 
d(P/PJ/de,  d2(P/Poo)/d0)  etc.,  from  the  tangent  wedge  approximation.  One 
approach  for  obtaining  the  cone  pressures  is  to  utilize  the  Mangier  trans¬ 
formation, 


d  M2  -Jc 

d8*/dx  =  ^|==-  2.20 

and  apply  this  to  the  equations  given  for  the  flat  plate.  An  alternate 
procedure  is  to  evaluate  the  coefficients,  d(P/Pso)/d6,  etc.,  from  the 
tangent  cone  approximation  and  provide  a  solution  directly  applicable  to 
the  cone.  This  has  been  done  for  K  >  1  so  that  the  first  coefficient  of 
d5*/dx  in  Equation  2.17  is: 
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and  a  considerably  longer  term  for  the  coefficient  of  the  second-order 
term. 


At  this  point,  an  analytic  iteration  procedure  Is  begun  by 

taking 


d  6  * 
dx 

equal  to  the  value  given  by  the  zero  pressure  gradient  boundary  layer 
solution  for  the  flow  over  a  cone.  This  Is  equivalent  to  using  the  flat 
plate  solution  with  the  application  of  the  Mangier  transformation  as  the 
zero  order  approximation  to 

d  6  * 
dx 

It  is  well  known  that  the  Mangier  transformation  is  applicable  only  when 
the  boundary  thickness  is  small  compared  to  the  body  radius.  Since  this 
restriction  is  violated  near  the  tip  of  the  cone,  even  in  absence  of 
interaction  effects,  the  validity  of  the  Mangier  transformation  under 
these  conditions  must  be  re-examined  in  the  future. 

Following  Lees  and  Probstein,  Equations  2.17  to  2.21  may  be 
used  to  give, 

f~  “  1  +  Fc  Xc  +  0  (xQ2)  2.22 

c 
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where, 


1  P» 


C  n/T  Mc  Pc 


2.23 


0=6. 


The  subscript,  0,  previously  used,  has  been  changed  here  to  c.  The 
value,  d,  is  still  given  by  2.2.  The  value,  Fc ,  is  plotted  in  Figure  2-6 
as  a  coefficient,  Fc  Mc/  "ydc  Moo,  versus  K.  The  coefficient,  Jj  /a/T, 
is  also  presented  for  comparison. 


^comparison  of  the  theoretical  methods  is  presented  in  Figure  2-7. 
The  ratio,  jj22-,  has  been  evaluated  from  the  compressible  flow  equations  as: 
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Since  a  strong  interaction  solution  is  not  available  (other  than  modifying 
the  results  for  the  wedge  by  the  Mangier  transformation),  the  liberty  has 
been  taken  _to  extend  the  first-order  weak  interaction  solution  to  large 
values  of  Xc*  By  analogy  to  the  case  of  the  wedge,  it  is  expected  that 
the  first-order  weak  interaction  solution  follows  the  strong  interaction 
solution  with  a  reasonably  constant  numerical  difference  of  percentage, 
since  both  solutions  are  linear  in  X-  Assuming  that  the  cone  solution  is 
correct,  the  errors  in  the  flat  plate  solution  with  Mangier' s  transforma¬ 
tion  (for  the  weak  interaction  case)  are  as  follows: 
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The  strong  interaction  solution  superimposed  on  Figure  2-8  (by 
applying  Mangier 's  transformation  to  the  flat  plate,  i.e.,  P/Pc  varying 
linearly  with  Xc/^0  >  *s  higher  than  the  cone  solution,  except  for  K  <  1. 
The  effects  of  slip  at  large  X  values  have  been  ignored  also,  which  would 
tend  to  reduce  the  induced  pressures  at  the  larger  values  of  X . 


Accordingly,  the  solution  as  given  by  Equation  2.22  will  be 
assumed  valid  over  the  whole  range  of  Xc- 
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WEAK  INTERACTION  REGION 


FIGURE  2—8 .  COMPARISON  OF  INDUCED  PRESSURES  ON  A  CONE 
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CHAPTER  3 

BOUNDARY  LAYER  GROWTH  AND  SHOCK  SHAPE  ON  SHARP 
CONES  AT  ZERO  ANGLE  OF  ATTACK 


3 . 1  COMPARISON  OF  ANALYTICAL  METHODS  FOR  BOUNDARY  LAYER  GROWTH 

Before  proceeding  to  the  cone,  the  case  of  the  sharp  flat  plate 
will  be  examined.  The  methods  to  be  compared  for  calculating  the  laminar 
boundary  layer  are  twofold  in  number: 

1.  Compressible  laminar  boundary  layer  theory,  and 

2.  Viscous  interaction  theory. 

Both  have  been  extensively  examined.  Van  Driest  has  numerically  solved 
the  laminar  boundary  layer  equations  using  the  Crocco  variables.  This 
solution  is  also  applicable  to  the  case  of  a  wedge,  if  the  independent 
parameters  are  interpreted  in  terms  of  the  local  values  at  the  edge  of 
the  boundary  layer. 

The  interaction  solution  should  give  better  results  than  that 
of  the  compressible  laminar  boundary  layer  theory  since  it  includes  the 
effect  of  the  pressure  gradient  and  pressure  level  with  the  interaction 
solution  (within  the  limits  already  discussed  for  xQ<l).  Intuitively, 
the  compressible  laminar  boundary  layer  theory  could  be  used  to  generate 
answers  identical  to  the  interaction  solution  by  considering  the  classical 
concept  of  the  displacement  thickness.  This  gives  a  new  equivalent  body 
from  which  to  calculate  a  new  shock  wave  and  resultant  property  changes 
across  the  shock.  An  interaction  procedure  can  be  devised  to  match  the 
boundary  layer  growth  (as  a  function  of  local  flow  properties)  to  the 
change  in  flow  properties  across  the  shock  wave.  The  interaction  theory 
offers  the  advantage  of  being  able  to  obtain  the  boundary  layer  thickness 
from  the  original  inviscid  flow  values  without  any  iteration  procedure 
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being  needed.  When  Che  boundary  layer  thickness  is  still  relatively  snail, 
however,  the  required  iterations  will  also  be  small,  since  the  shock  dis¬ 
placement  is  small.  The  principle  deviation  of  the  results  of  item  1  from 
item  2  of  Section  3.1,  would  then  be  expected  to  be  the  Increase  in  local 
Reynolds  number  due  to  the  induced  pressures  and  any  pressure  gradient  ef¬ 
fects  (Reference  5).  The  boundary  layer  thickness  may  then  be  obtained 
approximately  from  Reference  7  by  replacing  Re  with  Re  P/Pq,  where  P/Pg 
is  the  induced  pressure  caused  by  viscous  interaction. 


A  case  of  practical  interest  for  comparing  methods  is  a  hypersonic 
re-entry  condition;  accordingly,  a  comparison  of  these  methods  for  estimating 
the  boundary  layer  growth  is  given  in  Figure  3-1  at  H»  *  20.  In  this  and 
all  subsequent  altitude  calculations,  the  1959  ARDC  standard  atmosphere  is 
used.  It  is  seen  that  the  agreement  between  the  interaction  theory  and  the 
results  of  Van  Driest  with  a  pressure  correction  is  very  good.  The  difference 
between  the  two  solutions  in  the  weak  interaction  region  is  due  to  the  use 
of  different  viscosity  temperature  relationships.  The  effect  of  Prandtl 
number  is  also  noted  in  Figure  3-1.  Reference  6  presents  solutions  for  a 
continuous  variation  of  boundary  layer  thickness,  pressure,  skin  friction, 
and  heat  transfer  from  the  weak  interaction  region  to  the  strong  interaction 
region.  The  broken  line  of  Figure  3-1  represents  the  limiting  case  of 
strong  interaction  as  obtained  from  Reference  6. 


Figure  3-1  indicates  that  the  effect  of  the  pressure  gradient  on 
boundary  layer  growth  is  small  since  the  Van  Driest  solution  is  for  zero 
pressure  gradient,  and  the  strong  interaction  solution  is  based  on  a  value 
of  the  pressure  gradient  parameter ,  ji  =  0.286.  This  observation  has  been 
previously  noticed  for  a  highly  cooled  body^ ,  when  the  gas  near  the  body 
surface  is  much  cooler  and  more  dense  than  that  in  the  rest  of  the  boundary 
layer.  The  strearawise  pressure  gradients  are  then  not  very  effective  in 
perturbing  the  velocity  and  enthalpy  profiles  from  the  no  pressure  gradient 
case.  The  conventional  compressible  laminar  boundary  layer  solution  should 
then  be  adequate,  taking  into  account  only  local  pressure  effects  on  local 
Reynolds  number;  and  in  agreement  with  previously  stated  expectations.  The 
close  agreement  at  the  high  altitudes  where  the  boundary  layer  becomes  very 
thick  also  indicates  that  the  range  of  application  of  the  boundary  layer 
equations,  which  were  derived  for  a  thin  boundary  layer,  is  much  wider  than 
usually  supposed.  The  close  agreement  must  otherwise  be  held  to  be  fortuitous 
or  else  a  cancellation  of  the  various  effects  must  occur. 


The  study  of  the  flat  plate  provides  a  foundation  for  studying 
more  general  bodies,  such  as  the  sharp  cone.  The  general  conclusions  re¬ 
garding  the  interacting  laminar  boundary  layer  on  the  flat  plate  are  also 
expected  to  apply  to  the  case  of  the  cone.  The  actual  use  of  flat  plate 
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calculations  modified  by  a  constant  such  as  obtained  from  Mangler'B  trans¬ 
formation  must  be  made  with  caution.  The  Mangier  transformation  is  one  that 
reduces  the  boundary  layer  equation  for  axisymmetric  flow  to  those  for  two- 
dimensional  flow.  An  example  is  the  conversion  of  the  integral  momentum 
boundary  layer  equation  over  a  cone,  to  that  of  the  flow  over  a  wedge.  If 
the  coordinate  system  for  both  the  two-dimensional  and  axisymmetric  cases 
is  taken  so  that  the  direction  of  one  dimension  is  parallel  to  the  surface 
of  the  body  (e.g.,  along  the  generator  of  a  cone),  and  the  other  perpendicu- 
lar  to  the  surface,  then  the  forms  of  the  momentum  and  energy  equations  for 
the  axisymmetric  flow  are  Identical  to  those  for  the  two-dimensional  case. 
Therefore,  the  Mangier  transformation  must  be  able  to  change  the  form  of 
the  continuity  equation  from  the  axisymmetric  to  the  two-dimensional  form 
without  affecting  the  form  of  either  the  momentum  or  energy  equation.  This 
has  previously  been  demonstrated^. 

The  transformation  equations  are: 


Where  L  «  constant  reference  length 

r  ■  local  radius  of  revolution 


The  barred  variables  refer  to  the  two-dimensional  flow  case  and  the  un¬ 
barred  variables  refer  to  the  axisymmetric  case.  The  boundary  layer  growth 
over  a  flat  plate  in  compressible  laminar  flow  without  viscous  interaction 
is  given  by, 
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Peue 


d  y 


3.3 
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For  a  cone.  Equation  3.1  gives, 

x  -  r2  x/3L2  3.4 

If  3.2  and  3.4  are  substituted  in  3.3, 

5 


Equation  3.5  is  identical  in  form  to  3.3  with  all  barred  variables 
eliminated.  Equation  3.5  means  that  the  flat  plate  equation  may  be 
applied  to  the  case  of  the  cone,  providing  the  same  local  properties  are 
used  and  the  results  are  multiplied  by  a  factor  1 / n/~ 3 . 

The  boundary  layer  equations  assume  a  priori,  however,  that  the 
boundary  layer  thickness  is  thin.  This  means  that  the  Mangier  transforma¬ 
tion  (in  the  sense  of  correcting  flat  plate  results  by  a  constant)  is  not 
valid  in  the  case  of  a  thick  boundary  layer,  since  the  additional  terms 
(in  the  boundary  layer  equations)  required  to  describe  the  flow  will  not 
allow  the  form  of  the  axisymmetric  equations  to  be  reduced  to  those  of 
the  form  of  the  two-dimensional  flow.  In  more  general  terms,  the  purpose 
of  the  Mangier  transformation  is  to  transform  one  set  of  differential 
equations  to  another  set  for  which  solutions  are  known,  which  is  true  of 
any  transformation.  It  is  thus  seen  that  Mangler's  transformation  would 
not  invalidate  equations  of  motion  for  a  thick  boundary  layer,  but  rather 
that  a  different  set  of  differential  equations  must  be  solved. 

If  the  effects  of  viscous  interaction  are  considered,  it  is  con¬ 
cluded  that  Mangler's  transformation  gives  results  identical  to  equation  3.5 
for  the  case  of  weak  interaction,  since  6a nr'1,  and  is 
thin.  In  the  strong  interaction  region,  however  6a 
of  Equations  3.1  and  3.2  does  not  reduce  the  expression  for  6  to  a  form 
identical  to  Equation  3.3.  This  is  expected,  since  the  strong  interaction 
region  is  associated  with  a  thick  boundary  layer.  It  is  thus  apparent 
that  in  the  case  of  an  axisymmetric  flow  over  a  cone  with  a  thick  boundary 
layer,  and,  in  particular,  in  the  strong  interaction  region,  the  Mangier 
transformation  does  not  reduce  the  form  of  the  axisymmetric  equations  to 
the  two-dimensional  form.  Furthermore,  it  can  be  shown  that  any  similar 


still  relatively 
x3'  ,  and  the  use 
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transformation  of  the  type  x  ■  — ~  I  r”  d  x,  n  f1  2,  does  not  transform 

L  J 

the  axisymmetric  continuity  equation  into  two-dimensional  form  (e.g., 
Reference  14).  It  remains  to  be  shown,  however,  what  the  error  is  in 
assuming  a  correction  factor  of  the  type,  l/*/3,  to  the  flat  plate  calcu¬ 
lation  of  the  boundary  layer  thickness  in  the  strong  interaction  region. 


In  order  to  gain  insight  as  to  what  an  appropriate  engineering 
method  would  be  for  computing  the  thick  boundary  layer  on  a  cone,  three 
different  methods  will  now  be  examined  and  compared.  These  methods  are 

1.  Extrapolation  of  solution  derived  for  thin  boundary  layer 
to  the  case  of  thick  boundary  layer  (e.g.,  Reference  7 
with  Mangier1 s  transformation,  or  Reference  15  ata  •  0). 

2.  Viscous  interaction  theory  (Equations  2.10  -  2.12) 

corrected  by  factor,  3  Re  -  Re_„. 

3  >  c  FP 

3.  Solution  of  boundary  layer  equations  including  terms  of 
transverse  curvature  (Reference  16) . 

Solution  1: 


The  flat  plate  results  of  Van  Driest  may  be  used  by  replacing 
At Te  |  aT 3  Re  and  integrating  to  find  the  local  properties.  As  in  the 
case  of  the  flat  plate,  an  approximation  may  be  made  to  eliminate  the  re¬ 


quired  iteration  procedure  by  taking  ~  >»/"Re  as  ^  ■/3_ Re~ P/F, ,  and  taking 

the  local  flow  properties  as  the  original  inviscid  values,  assuming  that 
an  independent  determination  of  the  induced  pressures  can  be  made.  Inter¬ 
action  theory  for  the  induced  pressures  with  the  Mangier  transformation  or 
the  first-order  weak  interaction  on  a  cone  should  be  applicable  since  the 
transverse  curvature  effect  has  only  a  weak  effect  on  the  induced  pressures 
(p.  364,  Reference  1).  An  alternate  solution  is  given  in  Reference  15, 
where  the  results  ata.  «  0  may  be  used.  These  results  for  the  compressible 
laminar  boundary  layer  growth  over  a  cone  are  similar  to  those  of  Van  Driest 
except  that  Mangier' s  transformation  was  utilized  in  the  derivation,  and 
a  linear  viscosity-temperature  relationship  was  used.  The  boundary  layer 
thickness  may  be  presented  as  : 
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The  right  hand  side  of  Equation  3.6  integrates  to  the  identical  value  of 
the  right  side  of  Equation  2.10,  with  TQrig  replaced  by  the  local  value  Te. 

This  result  is  to  be  expected  since  the  assumptions  are  the  same  for  both 
derivations.  These  involve  the  flow  being  described  by  the  Blasius 
equation,  with  the  Mangier  transformation  utilized  in  the  case  of  the  cone. 
This  solution  also  represents  the  zero  order  term  for  the  case  where 
transverse  curvature  is  included,  and  will  be  examined  further  under 
Solution  3. 

The  problem  of  what  procedure  to  follow  in  the  iterated  solution 
will  be  examined  under  Solution  3.  The  procedure  for  Solution  1,  however, 
is  arbitrarily  taken  as  the  tangent  cone  method  using  NACA  Report  Number 
1133  to  obtain  the  shock  angle  (or  local  shock  slope) . 

Solution  2: 


For  the  viscous  interaction  solution,  modification  of  Equations 
2.10-2.12  gives  (Pr  -  1): 


a)  Weak  Interaction 


6*-W  3 Re 

x  V  C 


T 

1.730  ~ 
c 


+  0.664 


3.7 


which  is  identical  to  Equation  3.6. 
b)  Strong  Interaction 


x 


-  f. 


1/2 


6  <Hw/Ht>  X 


/  31/4  M 


3.8 


Solution  3: 


This  solution  is  taken  to  be  the  reference  solution,  since  it 
contains  all  the  aspects  of  the  actual  three-dimensional  flow,  within  the 
applicability  of  the  boundary  layer  approximations.  The  nature  of  the 
problem  is  such,  however,  that  an  iterative  procedure  must  again  be  devised 
to  match  the  boundary  layer  growth  with  the  inviscid  flow  properties  occurring 
across  the  actual  shock  wave.  The  examination  of  the  compressible  laminar 
boundary  layer  then  gives  a  solution  applicable  throughout  all  regions  of 
viscous  interaction  (i.e.,  weak  and  strong),  including  a  smooth  transition 
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between  the  two  regions,  and  accounting  for  the  axially  symmetric  boundary 
layer  growing  not  only  in  thickness  with  distance  along  the  surface,  but 
also  spreading  circumferentially  as  it  prows.  The  applicable  axially 
symmetric  boundary  layer  equations  are1®: 


Continuity , 


dy)  ur) 


&  fvr) 


3.9 


Momentum, 


dP 

dx 


dr  du 

dy  dy 


3.10 


Energy , 


_L 

Pr 


r 


dr  dh 


3.11 


Equations  3.9  to  3.11  were  transformed  in  Reference  16  to 
almost  two-dimensional  form,  and  an  explicit  expression  for  the  boundary 
layer  growth  derived  for  Pr  *  1 ,  neglecting  the  pressure  gradient.  The 
neglect  of  the  pressure  gradient  in  the  solution  of  the  boundary  layer 
equations  is  not  considered  to  be  a  serious  limitation  to  the  solution  if 
the  proper  pressure  level  is  accounted  for  (reference  flat  plate  results). 
This  is  taken  into  account  in  the  iteration  procedure.  In  addition,  the 
solution  accounting  for  transverse  curvature  is  presented  only  to  first 
order,  and  Probstein  and  Elliott  emphasize  that  the  first-order  effect  of 
transverse  curvature  and  a  self-induced  type  of  pressure  gradient  are 
independent  and  can  be  calculated  separately. 

To  first  order,  the  solution  is  then  as  follows: 

&*/£>*  =  1  +  '  (F./F  -  F  /2)  + . ,  3.12 

M  i  o  o 


where  5*  is  the  Mangier  value  representing  the  zero-order  solution  as 
obtained  from  Equation  3.7,  and  the  quantity  £(Fj/F  -  F  /?..)  is  the 
first-order  correction  accounting  for  the  effects  of  the  transverse 
curvature.  This  correction  is  presented  in  Reference  16  as  a  series 
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of  Integral  expressions,  which,  in  turn,  are  expressed  in  terms  of  the 
zero-order  and  first-order  boundary  layer  profile  of  velocity  and 
temperature.  By  specifying  a  specific  geometry  (viz.,  a  cone)  the  first- 
order  correction  has  been  evaluated  by  planimeter  using  graphical  inte¬ 
grations.  To  the  accuracy  of  this  type  solution,  the  results  are  as  follows: 

T 

F0  -  1.721  ^  +  0.664  Me2  3.13 

e 


T  T 

F  -  (0.330  +  1.277  +  0.218  mJ2)  ^ 

e  e 


+  (0.517  +  0.913  =r  +  0.233  M  2  )  -  M  2  3.14 

T  z  e  i  e 

e 


i 


3Re  tan  o 


3.15 


The  correction  has  been  calculated  for  y  =  1.4  (Pr  »  1)  and  given  in 
Figure  3-2  and  Figure  3-3  as  (&*/&*  -  1)  1/f,  versus  M  for  values  of 

r  - 1  2 

T  /T  «  0,  1.0,  and  (1  —  M  ).  Heat  transfer  plays  a  very  important 

we  i  e 

role,  as  well  as  Mach  number,  in  determining  the  correction  factor. 

Probstein  and  Elliott  examined  the  case  of  the  insulated  wall  at  M  *  3.16. 

g 

They  found  a  value  of  +0.073  for  the  ordinate  in  Figure  3-3  (compared  to 
0.03  in  the  present  calculations,  which  is  considered  fair  accuracy  for 
the  calculations)  in  contrast  to  what  might  be  expected  from  a  favorable 
pressure  gradient  analogy.  They  did  observe  that  when  the  effect  of 
compressibility  is  considered,  the  correction  does  become  negative  for  a 
sufficiently  large  Mach  number.  This  agrees  with  the  present  calculations. 
They  apparently  failed  to  observe,  however,  the  important  effect  of  heat 
transfer.  For  Tw/TAw  <  <  1,  it  is  possible  for  the  correction  to  increase 

with  increasing  local  Mach  number.  Physically,  this  may  be  explained  by 
the  very  high  densities  near  the  wall  in  the  case  of  the  highly  cooled  body. 
The  cold  wall  distorts  the  velocity  profile  to  increase  the  boundary  layer 
thickness,  insofar  as  the  effects  of  transverse  curvature  are  concerned. 

The  Mangier  thickness,  5*  ,  however,  still  exhibits  a  characteristic 
decrease  as  the  wall  is  cooled. 
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FIGURE  3-2.  TRANSVERSE  CURVATURE  CORRECTION  FOR  BOUNDARY  LAYER  THICKNESS 
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Concerning  the  iteration  procedure,  the  question  arises  as 
to  exactly  what  value  of  the  equivalent  body  should  be  used.  The 
usual  procedure  is  to  use  the  tangent  cone  method,  where  an  equivalent 
cone  is  formed  locally  at  each  axial  position  by  increasing  the  actual 
cone  angle,  a  ,  by  the  amount,  tan’^ddtydx) .  The  slope  of  the  dis¬ 
placement  thickness  is  thus  required  to  find  each  successive  iterated 
body.  The  viscous  interaction  effects  are  accounted  for  by  using  the 
local  flow  properties  in  the  calculation  of  the  slope  of  the  displace¬ 
ment  thickness.  As  in  Solution  1,  the  tangent  cone  method  will  also 
be  followed  here.  The  local  slope  of  the  displacement  thickness  will 
be  obtained  by  differentiating  the  applicable  equation  for  the  displace¬ 
ment  thickness;  then 


d5* 

dx 


2  (6VX) 


3.16 


n  equiv.  body  =  tan 


-1 


do*  ! 


dx 


+  0 


3.17 


The  determination  of  the  shock  slope  and  standoff  distance ,  however, 
will  be  obtained  in  this  solution  by  using  first-order  hypersonic  small 
disturbance  theory,  instead  of  NACA  Report  Number  1135.  The  geometry 
of  the  flow  regions  can  then  be  described;  as  well  as  allowing  a 
comparison  of  iteration  procedures.  The  comparison  will  be  a  direct 
one,  inasmuch  as  all  solutions  are  essentially  the  same  in  the  weak 
iteration  region.  The  equations  and  procedures  for  the  hypersonic  small 
disturbance  theory  will  be  presented  later. 


For  purposes  of  comparison,  a  case  of  practical  interest  was 
again  chosen:  0  *  10  degrees,  1^,  =  20,  Tw/T^  *  1.0,  as  a  function  of 

altitude.  The  following  designations  will  be  used  to  identify  the 
results : 

Solution  la  -  Van  Driest  with  pressure  correction  (Lees  and  Probstein 
first-order  weak  interaction  solution) 

lb  -  iterated  Van  Driest 

2a  -  weak  interaction 


2b  -  strong  interaction 
3  -  iterated  reference  solution 
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The  results  of  the  calculations  for  the  boundary  layer  growth  are  presented 
in  Figures  3-4  and  3-5.  Real  gas  effects  were  also  considered  in  Solutions  lb 
and  3  in  evaluation  of  the  property  changes  across  the  shock  wave. 

As  in  the  case  of  the  flat  plate,  the  iterated  and  pressure  corrected 
compressible  laminar  boundary  layer  solutions  of  Van  Driest  are  in  close  agree¬ 
ment  with  the  weak  interaction  theory,  within  the  region  of  its  applicability. 
The  reference  solution  is  somewhat  lower  in  this  region,  being  very  close 
at  150  kilofeet,  but  about  11  and  12  percent  lower  than  solutions  la  and  lb 
at  200  and  250  kilofeet,  respectively.  This  difference  is  apparently  due, 
for  the  most  part,  to  the  viscosity  value  assumed  in  the  boundary  layer 
equations.  A  slight  difference  is  also  Accounted  for  by  the  method  for 
determining  the  local  shock  slope,  The  difference  between  the  solutions  la 
and  lb,  and  the  reference  method  would  have  been  reduced  to  9.4  and  9.3  per¬ 
cent  at  200  and  250  kilofeet  respectively,  if  NACA  Report  Number  1135  had 
been  used  instead  of  hypersonic  small  disturbance  theory.  The  difference 
between  the  reference  method  and  interaction  theory  is  due  to  the  iteration 
procedure.  The  boundary  layer  thickness  would  be  reduced  by  the  iteration 
if  this  solution  is  considered  to  be  a  first  approximation  to  the  case  of  a 
thickened  boundary  flow  due  to  viscous  interaction.  This  reduction  occurs 
because  of  the  favorable  (negative)  pressure  gradient.  Order  of  magnitude 
estimates  of  the  reduction  based  on  Equation  2.13,  (see  p.  348  of  Reference  1)  , 
are  in  close  agreement  with  Figure  3-4. 

As  the  viscous  interaction  effects  increase,  the  weak  interaction 
solution  diverges  strongly  from  the  other  solutions.  Solutions  la,  lb  and 
3  agree  reasonably  well  in  the  range  studied,  but  the  strong  interaction 
solution  (2b)  is  low.  The  matching  point  for  the  interaction  solution  occurs 
between  250  and  300  kilofeet  for  the  10  degree  cone  at  x  ■  1  foot.  The 
effect  of  the  transverse  curvature  would  be  expected  to  be  felt  at  this 
point.  The  results  of  Figure  3-4  tend  to  verify  this  statement.  It  is  noted 
that  the  calculations  show  the  effects  of  transverse  curvature  to  be  negli¬ 
gible  in  the  weak  interaction  region,  and  vigorously  manifested  in  the  strong 
interaction  region  soley  at  350  kilofeet.  Solution  2b  is  lower  than  the 
reference  method  at  350  kilofeet  by  about  46  percent.  It  appears  that  in 
the  case  of  slender  cones  with  the  thickened  boundary  layers  occurring  with 
large  viscous  interaction,  Solution  2b  is  not  valid.  It  is  also  noted  that 
the  transverse  curvature  does  not  alter  the  local  body  slope,  but  causes 
only  a  uniform  change  in  thickness;  therefore,  it  thus  produces  no  change 
in  the  induced  pressures. 

The  same  general  comments  also  apply  to  the  comparison  at  values 
of  x  ft  1,  as  seen  in  Figure  3-5. 
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FIGURE  3-4.  COMPARISON  OF  VARIOUS  METHODS  FOR  COMPUTING  DISPLACEMENT  THICKNESS  ON  A  CONE 
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Another  indication  of  the  validity  of  the  methods  compared  for 
the  boundary  layer  thickness  is  to  compare  the  induced  pressures  predicted 
by  the  iteration  solutions.  These  results  are  presented  in  Figures  3-6 
and  3-7.  Away  from  the  vertex,  all  solutions  agree  with  the  first-order 
interaction  theory  of  Lees  and  Probstein  to  300  kilofeet.  At  350  kilofeet, 
solution  lb  falls  considerably  below  the  other  two  methods  at  x  M  foot. 

Near  the  vertex,  the  agreement  is  not  as  good  (Figure  3-7).  The  main 
difference  between  solutions  lb  and  3  is  due  to  the  equivalent  body  slope 
as  obtained  from  the  boundary  layer  thickness  calculations.  This,  then, 
can  be  traced  back  to  the  assumption  made  for  the  viscosity-temperature 
relationship.  The  accuracy  of  any  of  the  methods  very  near  the  vertex, 
however,  or  at  350  kilofeet,  is  open  to  considerable  question;  due  to  the 
neglect  of  other  important  considerations,  aich  as  slip  flow,  nonequilibrium 
flow  effects,  etc.  The  calculations  at  these  conditions  are  primarily  of 
academic  interest,  and  to  indicate  general  trends. 

Some  experimental  results  are  available  in  Reference  17  for  in¬ 
duced  pressures  on  3°  semi-vertex  angle  cones.  Sample  calculations  by 
Solution  3  have  been  carried  out  to  compare  with  this  data  at  one  condition. 
These  results  are  given  in  Figure  3-8.  Since  Mw  0|,  <  1,  the  interaction 
theory  is  not  valid  for  these  test  data.  The  test  models  were  cooled  to 
thermal  equilibrium  and  then  the  pressure  measurements  were  taken.  The 
models  were  therefor,  at  or,  slightly  greater  than  the  adiabatic  wall  tempera¬ 
ture.  The  induced  pressure  increment  (P2  -  Pc)/Pc  is  very  low.  It  was  found 
that  the  magnitude  of  the  increment  at  these  low  levels  was  extremely  sensi¬ 
tive  to  the  accuracy  of  the  calculations.  Since  the  calculations  were  carried 
out  by  hand,  the  deviation  of  prediction  and  data  is  within  the  possible 
error  involved  in  the  calculations. 

On  the  basis  of  the  comparison  just  made,  and  from  the  standpoint 
of  ease  of  calculation,  it  is  recommended  that  the  boundary  layer  growth  on 
a  cone  be  computed  as  follows: 

1.  Solution  2a  (weak  interaction)  may  be  used  in  the  weak 
interaction  region. 

2.  Solution  la  (Van  Driest  with  Lees  and  Probstein 
first-order  weak  interaction  pressure  correction) 
may  be  used  to  extend  solution  2a  into  the  strong 
interaction  region. 

The  first-order  weak  interaction  solution  of  Lees  and  Probstein 
for  predicting  the  induced  pressures  should  be  appropriate.  It  is  felt  that 
these  results  for  the  boundary  layer  growth  and  induced_pressures  should  be 
applicable  to  cones  with  a  ^  10  degrees  in  the  range  5  <  a  <  20  degrees. 
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FIGURE  3-6.  COMPARISON  OF  INDUCED  PRESSURES  AS  A  FUNCTION  OF  ALTITUDE 
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FIGURE  3-7.  COMPARISON  OF  PRESSURE  DISTRIBUTIONS  ON  A  CONE 
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FIGURE  3-8.  INDUCED  PRESSURE  INCREMENT  ON  A  CONE  (DATA  OF  REFERENCE  17) 
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3.2  DETERMINATION  OF  EQUIVALENT  POWER  LAW  BODY  AND  SHOCK  SHAPE. 

The  growth  o£  the  boundary  layer  on  a  cone  has  already  been  de¬ 
termined  as  6*(x‘),  where 


Ox,n 


o 


We  may  now  show  that  a  boundary  layer  growth  of  the  type  described  over 
this  sharp,  slender  cone  produces  a  new  equivalent  body  characterized  by: 


y 


equiv 


r.O!xm,  l/2'rm<l 

D  — 
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The  rotation  of  the  (x1 ,  y’)  coordinates  into  the  (x,  y)  coordinate  system 
is  given  by: 


x 

y 


X*  COS  0 

x1  sin  o 


sin  o 
cos  CJ 
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In  the  (x1,  y’)  coordinate  system,  the  boundary  layer  growth  is: 


6*=  y’ 


K  x 


.n 


1/2  <  n  <  3/4 


3.20 


It  can  be  shown  that  0  <  K  <  1,  where  K  _»0  as  Re  -» °°.  In  the  (x,  y)  co¬ 
ordinate  system,  the  substitution  of  Equation  3.18  and  3.13  gives: 


v  .  ■  r. 

equiv  b 


Kx 


n 


1-n 
cos  o 


1  +  —  tan  a 


+  tan  o 


3.21 


3-20 


9Sr^J^g(ofSorn/umy-, 

AERONUTRONIC  DIVISION 


For  slender  bodies  with  tan  o<<1,  n  <  1, 


Kx 
1-n 
cos  a 


1  +  n  —  tan  o  J  +  x  tan  o 
x 


3.22 


and 


1-n  2-n  1 -q 

m  cos  o  tan  ox  +  K  x  m  T(x  cos  o )  tan  o  +  K]  x 

b  1-n  1-n  ,  _ . 1-n  _ 

_  „  (x  cos  a )  -  K  n  tan  a 

cos  o  x  -Kn  tan  a 


3.  23 


It  Is  now  postulated  that: 


rb  a  x 


3.24 


An  empirical  correlation  was  made  to  determine  the  constant  of  proportionality 
in  3.24  by  comparison  to  3.23.  This  results  in: 

0.30 

rb  (tan  cr  +  1.10  K)  xl*<x-n>K  3  25 


The  percent  deviation  of  the  empirical  value  of  rj,  from  that  of  the  exact 
value  of  rb  is  given  in  Figure  3-9  for  several  values  of  the  determining 
parameters.  The  deviation  ranges  from  -13  to  +207.  at  x  =*  0.1. 

The  accuracy  of  a  point  so  near  the  leading  edge  is  somewhat 
questionable.  The  region  of  primary  interest  is  considered  to  be  0.5  <  x  <  2, 
where  the  agreement  is  approximately  +  5  percent. 

The  equations  of  motion  for  inviscid  hypersonic  flow  over  slender 
bodies  can  be  reduced  to  simpler  form  by  incorporating  the  hypersonic-slender- 
body-approximat ions .  This  use  of  small  disturbance  theory  in  hypersonic  flow 
has  been  summarized  by  Van  Dykel8.  The  reduced  equations  are  valid  provided 
6  «1  and  1/M  oo6<0(l),  where  6  is  a  characteristic  body  or  streamline  slope, 
and  Moo  is  the  freestream  Mach  number.  Approximate  analytical  solutions  for 
the  inviscid  hypersonic  flow  over  power  law  bodies  of  the  type  r(,  Ct  x"1  has 
been  given  by  Mirels^^.  The  resulting  equations  give  the  zero-order  and 
first-order  hypersonic  flow  over  those  blunt-nosed  configurations  whose 
zero-order  shock  shape  follows  a  power  law  variation.  If  the  classical  con¬ 
cept  of  the  Prandtl  displacement  thickness  is  considered,  it  is  seen  that  in 
the  high  altitude  high  Mach  number  region  of  interest,  the  thickened  boundary 
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FIGURE  3-9.  DEVIATION  OF  EMPIRICAL  BODY  RADIUS  FROM  EXACT  SOLUTION 
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layer  growth  on  a  sharp  cone  produces  an  equivalent  body  which  Is  effectively 
blunted.  Since  the  boundary  layer  growth  can  be  expressed  as  a  power  law 
variation  of  the  cone's  slant  length,  It  Is  to  be  expected  that  the  equiva¬ 
lent  body  can  be  represented  as  a  power  law  body.  This  has  just  been  verified 
by  an  empirical  correlation.  Within  the  accuracy  of  the  hypersonic  small 
disturbance  theory,  a  method  Is  available  to  estimate  the  entire  lnvlscid 
flow  field  downstream  of  the  vertex,  assuming  of  course,  that  the  equivalent 
body  Is  known,  a  priori.  This  has  been  determined  to  be  the  case.  In  the 
following  discussion,  only  the  applicable  equations  and  instructions  for 
their  use  will  be  presented,  since  a  complete  description  is  adequately  given 
in  References  18  and  19. 


The  shock  shape  is  given  by: 


!R_  „  i_ 

1  i  _ 

W  ( 

,  2(l-m)‘ 
x  | 

rb  ’  % 

1  *r 

<M»v2  ' 

L  ' 

slope  is: 

dR  _  1 

drb 

- 

i  +  — 
m 

W  I , 

dx  ”  qb 

dx 

(Moo&b)2  1  1 

2(l-m) 


or 


dR 

_ 9 

dx 


dx 


1  1_  6* 
2  nb  x 
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3.27 


3.28 


The  quantities  %  and  a^  are  given  in  Figure  3-10  as  functions  of  V  and  the 
pressure  gradient  parameter,  0 .  The  characteristic  slope,  &b.  Is  for  the 
slender  bodies  under  consideration: 
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b 


x 

L 


tan 


5* 

a  +  —  /cos  o 


x»L 


-  tan  o  +  (5*/x)  T 
x“L 
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The  snalysis  of  References  1,  2  and  19,  indicates  that  for  many  considera¬ 
tions  the  variation  of  properties  normal  to  the  surface  may  be  negligible 
within  the  accuracy  of  the  engineering  approximations  employed.  As  an  example, 
the  constant  density  approximation  used  in  Reference  2  proved  highly  succes- 
ful,  and  the  variation  of  pressure  as  obtained  from  Reference  19  was  found  to 
be  negligible  within  approximately  5  percent  for  the  cases  of  Interest  where 
a  relatively  thin  shock  layer  (or  invlscid  region)  exists.  Accordingly,  the 
present  analysis  will  assume  no  variation  in  pressure  or  density  (consequently, 
temperature  as  well)  normal  to  the  surface  and  utilize  the  values  Immediately 
aft  of  the  shock  wave  in  the  invlscid  flow. 
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CHAPTER  4 

FLOW  FIELD  DETERMINATION  OVER  6.33 
AND  10  DEGREE  SEMI-VERTEX  ANGLE  CONES 


4.1  BOUNDARY  LAYER  AND  SHOCK  LAYER  PROFILES 

With  the  previous  results  available,  we  are  now  in  a  position 
to  estimate  the  geometry  of  arbitrary  conical  flow  fields.  The  temperature 
and  density  profiles  therein  can  be  obtained  directly  in  the  shock  layer, 
but  have  not  been  ascertained  for  the  boundary  layer.  A  simple  estimate 
of  the  boundary  layer  profiles  may  be  made  by  using  the  Crocco  integral 
expressed  in  terms  of  enthalpy^»2,4. 


where  velocity  is  considered  to  be  the  independent  variable.  This  result 
follows  directly  from  a  simultaneous  solution  of  the  momentum  and  energy 
equations  for  Pr  »  Le  ■  1,  dP/dx  ■  0,  and  irrespective  of  the  body  shape. 
With  a  constant  pressure  across  the  boundary  layer  the  temperature  and 
density  profiles  may  be  determined  from  Equation  4.1  if  the  boundary  con¬ 
ditions  are  known.  These  have  been  determined  above.  In  those  cases 
where  dissociation  become  important,  the  temperatures  may  be  obtained  from 
Reference  20  as  T  ■  §  (h,  P2). 

In  hypersonic  flow  without  blowing  the  velocity  profile  becomes 
nearly  linear^.  With  this  approximation  and  Equation  4.1,  a  simple  solu¬ 
tion  is  available  to  estimate  the  entire  flow  field.  The  assumptions  of 
zero  pressure  gradient  and  linear  velocity  profile  are  not  considered  to 
be  serious  limitations  to  the  accuracy  of  the  predicted  results.  As  a 
check  on  the  possible  error  involved  in  the  zero  pressure  gradient  as¬ 
sumption,  a  comparison  of  the  enthalpy  distributions  through  the  boundary 
layer  was  made  at  P  ■  0  and  P  ■  0.286.  The  latter  condition  corresponds 
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1  -1/2 
to  a  locally  hypersonic  flow  with  the  boundary  layer  growing  as  x 

This  comparison  is  given  in  Figure  4-1.  The  distribution  for|3-  0.286 

were  obtained  from  the  calculations  of  Li  and  Nagamatsu  presented  in 

Reference  21' for  a  similar  solution  of  the  boundary  layer  equations. 

Calculations  have  been  made  for  the  flow  fields  over  two  cones 
with  6.33  and  10  degree  aemi-vertex  angles.  The  computations  were  carried 
out  for  a  typical  re-entry  condition  of  ■  25,  at  an  altitude  of  250 
kilofeet.  These  results  are  presented  in  Figure  4-2  to  4-7. 


4.2  HEAT  TRANSFER  AND  SKIN  FRICTION  ON  CONES 

In  a  manner  similar  to  that  followed  for  boundary  layer  thick¬ 
ness,  the  skin  friction  and  heat  transfer  on  sharp  slender  cones  may  be 
evaluated.  It  was  shown  that  a  smooth  transition  from  the  weak  to  strong 
interaction  region  may  be  obtained  by  neglecting  pressure  gradient  effects 
and  considering  only  the  local  properties.  This  may  be  approximately 
accounted  for  by  considering  the  change  in  local  properties  due  to  the 
induced  pressures.  In  the  case  of  skin  friction  and  heat  transfer  however, 
the  transverse  curvature  effects  may  not  be  negligible! , 16  and  should  also 
be  included  in  a  general  analysis. 


Following  this  procedure  and  neglecting  slip  effects,  the  local 
skin  friction  coefficient  on  a  sharp  slender  cone  is  given  by: 

4.2 


Cf  -  0.664 

\ 

J  3C  P/Pc 

T 

1  +  £,  (0.517  +  0.913  y  +  0.121  (y-l)M2 

/  Re. 

L 

c 

With  P/Pc  given  by  Equation  2.22  the  mean  skin  friction  coefficient  is: 


C--  2~  *  ^■443Vi?  <1+  rcXc,3/2  + 
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2  RrSfe  (0.517  +  0.913  ■—  +  0.121  (7-1)  M2) 
L  c 
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FIGURE  4-1.  COMPARISON  OF  ENTHALPY  DISTRIBUTIONS  FOR 
EFFECT  OF  PRESSURE  GRADIENT 
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FIGURE  4-2  .  BOUNDARY  LAYER  THICKNESS 
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FIGURE  4-3. A.  SHOCK  LAYER  TEMPERATURES 
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3/2 

U  +  VO  Fc*c  1/2  <Fc*c> 

- -  -T*  (1  +  W  '  -  -V 


1/2 


In 


1/2' 

1  +  (1  +  F  X  ) 

(F  7  )1/2 
(Fc  *  c>  J  j 


4.3 


Equations  4.2  and  4.3  now  include  both  the  viscous  interaction 
and  transverse  curvature  effects,  without  considering  them  separately.  An 
alternate  procedure  for  estimating  the  skin  friction  forces  is  to  obtain 
a  value  by  summation  of  the  individual  effects.  With  the  substitution  of 
the  appropriate  boundary  conditions  the  local  skin  friction  is  obtained  as: 


C£-  (C£) 


Mang 


(Cf) 


Mang 


T  F  \ 

1.152  +  0.733  ~  +  (0. 5287-0. 126)M  2  -£"4  4.4 


As  mentioned  previously,  transverse  curvature  might  have  a 
significant  effect  on  the  local  skin  friction.  The  first-order  transverse 
curvature  correction  to  the  local  skin  friction  is  given  by  Probstein  and 
Elliot  in  Reference  16  as: 


cf  •  <cf> 


<c£> 


Mang 


Mang 


0.517  +  0.913  j1-  +  0. 121(y-l)M' 
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2j - 1 —  \p£ 

cjy3Tan0  V  (Rx) 
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The  total  local  skin  friction  including  interaction  and  transverse  curvature 
is  now: 


,  /  3C 

“  0.6WV7E 
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C  0  T 

(0.517  +  0-913  rr  +  0.12Uy-l)M*  X 

>/,  „3  Tc  c  c 


4.6 


r  L  c 

i 

♦H 

0.528  7  *  0.126  + 

[  < 

T 

w 

T 

c 

1.152  1 

<  J 

)  \ 

4-10 


AERONUTRONIC  DIVISION 


We  may  express  the  local  heat  transfer  rate  (by  utilizing 
Reynold's  analogy)  In  terms  of  a  Stanton  number: 


l  2 

sti  -  1  s  V? 


where , 


St,  5 


/  _ 

p* J  e  °P 
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4.8 


and  the  local  heat  transfer  rate  Is  given  by: 


St  p  u  (H  -  H  ) 
£  r  e  Aw  w 
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4 . 3  AXIAL  FORCE  COEFFICIENT  AT  ZERO  ANGLE  OF  ATTACK 

In  order  to  check  the  above  results,  a  comparison  with  experi- 
mentaljiata  for  axial  force  coefficients  measured  in  a  hot  shot  tunnel22 
at  =  20  has  been  made.  The  axial  force  coefficient  is  obtained  by 
adding  the  individual  contributions  for  the  axial  force: 


a-  0 


pF  B  rN  P 
°A  +  CA  +  CA  +  CA 
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The  force  components,  C^,  are  estimated  by: 


F  A 


cos  a  (integrated  skin  friction  drag) 
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(flat  nose  drag) 
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c 


p 

A 


1  m2 

2  ”00 


(1  +  3  7FCTC)>1 


(pressure  drag  from 
Equation  2.22) 


Results  for  this  comprison  is  presented  in  Figure  4-8.  The 
base  pressure,  Pg>  has  been  estimated  from  Reference  23.  The  contribu- 
tions,  C®  and  cjj  are  negligible  at  the  lower  Reynolds  numbers.  The  max¬ 
imum  deviation  is  about  +  10  percent.  It  is  felt  that  the  analytical 
estimates  are  adequate  within  the  accuracy  of  the  assumptions  involved. 

The  principal  assumptions  to  be  considered  are  the  neglect  of  slip  effects 
and  the  approach  employed  in  using  Equation  4.10.  It  is  also  noted  that 
is  very  sensitive  to  the  accuracy  of  the  calculations  and  the  methods 
employed  for  computing  the  inviscid  flow  values.  As  an  example,  the  use 
Re„ 

of  - —  ■  1.20  (instead  of  the  value  used,  e.g.,  1.50)  increases  C.  by 

Reoo  4  A 

about  14.3  percent  at  Reo^"  1  x  10  .  For  this  wind  tunnel  data,  the 

calculations  are  especially  sensitive  to  the  viscosity.  At  the  low  ternp- 
peratures  involved,  Sutherland's  law  is  inadequate,  and  data  was  used  from 
NAVORD  Report  1488  (Vol.  5). 

Comparison  with  experimental  data  has  also  been  made  by  using 
the  summation  method  for  obtaining  the  skin  friction  forces: 


C_  •  1.53  cot  6 
°F 


A  v/JL 

V(\) 


i  + 


Cote 


0.446  +  0. 788  “  +  0.105  (Y-1)M2^(XL)  + 
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3  F  . .  0.733  w  1.152.,-  . 

+  2  7  (0.528  7-0.126  +  j-  —  -  j  )(xL) 

M  c  M  c 

c  c 
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At  the  lower  Mach  numbers  (M<10)  the  base  drag  contribution  can 
be  significant  and  it  is  necessary  to  know  the  value  of  the  ratio  Pb/Fo,  • 
This  is  a  rather  complex  calculation  since  it  appears  to  be  a  function  of 
the  local  flow  conditions  at  the  base  of  the  vehicle  as  well  as  the  local 
boundary  layer  characteristics^. 
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FIGURE  4-8.  COMPARISON  OF  THEORY  AND  EXPERIEMENT  FOR  AXIAL  FORCE  COEFFICIENT 
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Before  concluding  this  discussion  of  Interaction  theory,  one 
final  point  of  interest  remains.  In  the  case  of  using  the  square  root 
of  the  ratio  of  the  interaction  pressure  to  the  inviscid  pressure  to  ob¬ 
tain  the  total  skin  friction  in  the  presence  of  boundary- layer- shock- 
wave-  in  terac  t  ion  : 


c 


where  the  subscript  i  refers  to  the  case  with  interaction. 

Bertram  has  shown  that  this  approach  gives  very  good  results 
in  the  case  of  a  flat  plate. 

Applying  this  procedure  to  the  case  of  a  cold  wall  cone  at  high 
Mach  numbers,  and  letting  Equation  2.22  represent  the  pressure  ratio: 


V1  *  P c  X  c 


4.12 


which  for  small  values  ofXc  is  approximately 


F 

c 


4.13 


For  these  same  conditions  Frobstein's  theory  gives: 


— -  1  +  (0.528 

Cf 


0.126 


F  X  • 
c  c 


4.14 


_  Assuming  aT»  1.4  there  is  approximately  a  12  percent  difference 
in  the  X  term  between  Equations  4.13  and  4.14. 


4-14 


AERONUTRONIC  DIVISION 


Experimental  drag  data  for  a  9°  cone  obtained  under  various  test 
conditions  are  shown  in  Figure  4-9.  Also  shown  in  this  figure  are  appropriate 
drag  curve  segments  calculated  by  means  of  the  Interaction  theory.  It  is 
apparent  that  reasonable  agreement  exists  between  the  calculated  curves  and 
the  test  points  at  the  higher  Reynolds  numbers.  However,  at  the  lowest 
Reynolds  numbers  tested,  the  disagreement  is  substantial.  With  the  limited 
test  data  available,  presented  in  this  form,  it  is  not  possible  to  establish 
precisely  where  the  theory  and  the  test  data  begin  to  diverge. 

Trends  similar  to  those  previously  stated  are  observed  for  a  6.33° 
cone,  as  shown  in  Figure  4-10. 

Interpretation  of  the  various  trends  shown  in  Figures  4-9  and  4-10 
will  be  left  to  a  later  section.  The  principal  point  to  be  made  here  is 
that  within  the  region  of  applicability,  the  interaction  theory  gives  re¬ 
liable  results. 

It  should  be  noted  that  the  NOL  data  presented  in  Figure  4-10  are 
of  a  reduced  order  of  accuracy.  These  data  were  obtained  in  a  ballistic 
range  facility  where,  due  to  the  low  density  conditions  under  which  the 
test  was  run,  it  was  not  possible  to  prevent  the  model  from  undergoing  a 
sustained  pitching  oscillation  as  it  traversed  the  range.  Therefore,  the 
actual  measured  drag  of  the  model  included  some  angle  of  attack  effects. 

The  data  presented  in  Figure  4-10represent  the  measured  data  reduced  to 
zero  angle  of  attack  conditions  by  means  of  an  extrapolation  technique. 

Because  of  the  inherent  uncertainties  of  this  procedure,  the  disagree¬ 
ment  between  the  NOL  data  and  the  theory  shown  in  Figure  4-9  is  not  con¬ 
sidered  significant. 

A  substantial  part  of  the  experimental  data  shown  in  Figures  4-9 
and  4-10 have  been  obtained  from  unpublished  sources.  Since  this  informa¬ 
tion  is  not  generally  available,  a  brief  summary  of  the  principal  conditions 
and  results  for  these  tests  is  given  in  Table  II. 
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FIGURE  4-9.  COMPARISON  BETWEEN  EXPERIMENT  AND  THEORY  9 
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FIGURE  4-10.  COMPARISON  BETWEEN  EXPERIMENT  AND  THEORY  6.3  CONE 
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CHAPTER  5 

APPROXIMATE  INTERACTION  THEORY  AND  EXPERIMENTAL  CORRELATION 


5.1  APPROXIMATE  INTERACTION  THEORY 

While  the  complete  interaction  theory  discussed  in  the  Section's  II 
III,  and  IV  indicates  agreement  with  experimental  data  it  is  somewhat  cumber¬ 
some  in  practice.  The  present  section  develops  an  approximate  form  of  the 
interaction  theory  which  permits  the  rapid  calculation  of  cone  drag  character 
istics  at  hypersonic  speeds.  This  simplified  form  of  the  interaction 
equation  also  provides  the  basis  for  a  correlating  procedure  which  unifies 
the  experimental  drag  data  into  a  single  curve. 

In  the  derivation  vrfnich  follows,  it  has  been  found  convenient  to 
express  the  flow  properties  in  terms  of  the  free  stream  Mach  number  and 
Reynolds  number  rather  than  the  conditions  behind  the  shock.  This  pro¬ 
cedure  introduces  certain  restrictions  to  be  noted  subsequently.  However, 
it  does  facilitate  the  direct  comparison  of  the  drag  of  various  cone  con¬ 
figurations  . 

The  leading  assumption  is  that  the  free  stream  Mach  number  is 

hypersonic,  i.e.,  M  >>1. 

00 

From  the  definition  of  X  it  follows  immediately  that: 


M 

For  convenience,  we  also  define,  F  =  F'  _ dc . 

C  M 

c 
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Assuming  the  linear  viscosity  temperature  law  and  a  slender  body,  the  following 
relations  apply: 


C 

t 

c" 


T 

1/4  M  \ 

00 

c  1 

T 

M 

c 

»  / 

1/2 


Me 

N. 


p 

_ £ 

P 


7/2 


Substituting  these  relations  Into  equation  5.1  we  obtain: 

\l/2 


M 

_ c 

M„ 


P 

OQ 

P 


5.2 


With  this  expression  for  X  ,  the  cone  pressure  drag  (including  vacuum  base 
drag)  becomes: 

«.  -  -v  f  »  ♦  4  r/4-r  -.*1  ,3 


TM 


For  large  Mach  numbers  K>>  1  and 


K  „  lirtiis 

s  r+3 


1  +  |  F’l 

fF  \1/2  I 

~~  d  X 

3  | 

[cl  C  “ 

Hence  _c  _  „  (r  +  1)  „2 

P  ~  8  r  ,  .  „,2  K 


(T  +  3) 


The  coefficient  F'  is 

F'  - 


•JT 

For  K  >>1,  it  can  be  shown  that: 
d(P/P  ) 


1  1  [  1  d(P/P~>  ' 

r  p/p»  Mcc  d  9 

«-  Je  =  9, 


_i_  _ 

m„  d  e 


e  =  ec 


16  r  111  l)  k 

(r  +  3)2 
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Therefore  F'  may  be  approximated  as: 


l 

K 


5.5 


Expressing  the  cone  angle  6  in  terms  of  the  body  length  to 


diameter  ratio 

(cot  e  -  2  l/D) 

and  assuming  a  slender  body,  Equations 

5.4  and  5.5  can  be 

rewritten: 

2 

P 

c 

*  2  y 

7+1 

M 

00 

p*>  1 

(7+3)2 

(  £/D)2 

F '  » 

4 

Jg/D 

V3 

M 

00 

Substituting  these  relations  into  Equation  5.3  for  the  cone  pressure 
drag,  we  obtain: 


CDP 


JLi  -1 . 

(  7+3)2 


5.8 


C 

00 


The  above  equation  contains  the  viscosity  proportionality  constant 
based  on  the  free  stream  conditions.  From  the  power  law. 


C 

00 


1/4 


The  ratio  Tref/^  may  be  evaluated  by  Eckert's  reference  en 
thalpy  method.  00 

Approximating: 


+ 


(1  +  M  )  by  M 

Z  oo  Z  oo 


and  expressing  the  body  wall  temperature  in  terms  of  the  stagnation 
temperature,  Eckert's  reference  enthalpy  method  gives: 


ref 


(Pr) 


1/2 


<¥>  «» 


,  0.5 
'(Pr) 


1/2 
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For  convenience  we  define: 


Thus : 


E 


<^)( 


0.5 

(Pr)1/2 


T 

~  +  0.22) 

A0 


c  - 

00 


*0 


1/4 


5.9 


5.10 


and  che  cone  pressure  drag  equacion  becomes: 


CDp 


7+1 

("Tfs)2 


(i/D)2  + 


d 

_ c_ 

7+3 


(eF/8 


5.11 


He  shall  leave  Che  pressure  drag  term  in  Chis  form  for  Che  Cime 
being  and  proceed  wich  Che  skin  friccion  drag.  IC  is  worCh  noCing,  however, 
ChaC  Che  firsc  Cerm  in  Che  above  equacion  is  Che  inviscid  pressure  contri- 
buCion,  while  Che  second  term  conCains  Che  inCeracCion  pressure  conCribuCion. 


The  skin  fricCion  drag  consisCs  of  three  terms,  considered 
separately;  the  Mangier,  interaction,  and  transverse  curvature  terms.  The 
first  term  is  the  Mangier  skin  friction,  the  second  represents  the  effect 
of  transverse  curvature,  while  the  third  represents  the  effect  of  the 
boundary- layer-shock-wave  interaction . 


Since , 


q  P 

OO  I 


M 


the  skin  friction  equation  may  be  written  in  the  form: 

_  -  3/4 
CDp  -  4.33  '/T’^75  /-v  1/8 


M 


+  6.12  ( L/D) ‘ 


+  10.59  (L/D)' 


all. 


5i' 

mc2JI<e>1/8  v^vtJ 

f  ^  If  M  3/4-1  2 

Ql  d  _ L  _(g_. 

L7  CJU1/8  V(VT J 
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where  0  -  (0.446  +  0.788  rr  +  0.105  ( 7-I)  M2) 


M 


®-  <0.528  7-0.126+  °-;7-33- 


M 


Combining  Equation  5.12  for  the  skin  friction  with  Equation  5.11  for  the 
pressure  drag  the  total  drag  Is 


■y+1  4 

- T2 


•SSL 


3/4 

i _ 

v«sr 


«>■  w  ”  **  d' + 

3/4 

t(W{612[f]+10'59[^^]}^fe- 


5.13 


We  now  proceed  to  simplify  the  terms  In  the  brackets.  The  term 
dc  may  be  rewritten,  under  the  assumption  of  large  Mach  number: 


dc  *  A(Pr) 


+  B  (Pr)  (7  -  1) 


Taking  A(Pr)  =  0.968  and  B(Pr)  -  0.145: 


d  =  (  7*  1)  (0.145  +  0.484  ) 

C  T0 


Thus : 


( 


d  +  4.33^7 

*ly  c 


J  =  5.61sT7  J 


w  0.226  0  754 

1  +  0.754  rr  - 

Lo  y  y 


Iwl 

ToJ 


It  is  now  proposed  that  the  last  two  terms  in  the  above  brackets 
be  numerically  evaluated  using  a  value  of  7  *  1.4.  This  gives: 


( 


=  4 . 7nT7  1  +  0.258  —  j 


d£  -  4 . 33n/"7  ^4.7nT7  |l  +  0.258 


5.14 


The  above  expression  is  correct  for  7*  1.4,  but  is  somewhat  In 
error  for  other  values  of  7  .  For  a  7  »  1.67,  and 
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there  is  approximately  a  3  percent  error  in  the  above  expression,  while  for 
the  conditions  ym  1.67  and 

T 

^  “  1 
T0 

the  error  is  approximately  10  percent. 

The  other  term  to  be  simplified  in  Equation  5.13  is: 


[-$-] 

+  10.59 

'  &  * 

l  M2  . 

c 

7 

) 


Expressing  the  body  wall  temperature  in  terms  of  the  stagnation  temperature: 


■Q-  0^6  +  0  394  (7_  x)  Ji  +  0.105  (7-  1) 

M  M  A0 

c  c 


At  large  Mach  numbers  the  first  term  is  negligible. 

Thus:  ^  (  y.  1}  (0  105  +  Q.394  -Si  ) 

0 


M 


©L 


Following  a  similar  approach  for  'aK-  d  : 

7  c 


5.15 


d  «  (7  -  1) 

7  C  ' 


0.0766  (1  -  +  0.30S  (1  -  )  + 

7  lQ  1 


0.177  1 

0 


(1 


i>] 


Substituting  the  value  7“  1.4  into  the  terms  within  the  brackets 
the  above  reduces  to: 

©  Tw  /  Tw  ‘  2 

^d  «  (  7  -  1)  0.0635  +  0.227  rr  +  0.0505  [—■ 

7  c  T0  l  Tq 


5.16 
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For  a 
»  3.5  percent. 


T 

~  -  1 
T0 


7  ■  1.67  and 
If  y  =  1.67, 


the  error  is 


T 

w 

—  ■  0  the  error  in  the  above  expression  is 
T0 


approximately  12  percent. 


Making 
5.16,  the  result 


use  of  the  approximate  relations  of  Equations  5.15  and 
is  obtained: 


10.59 
2  , 


7 


(  7-  D 


T 

1.313  +  4.81  r2  + 
T0 


5.17 


For  most  cases  of  practical  interest  (i.e.,  atmospheric  re-entry 
of  a  flight  vehicle)  the  temperature  ratio 


T 


w 


is  quite  small.  Even  for  an  adiabatic  wall 


2 

^0.72, 


and  since  the  last  term  in  the  brackets  above  is  only  8  percent  of  the 
previous  term  this  last  term  will  be  neglected. 

The  total  drag  for  the  cone  is  now  expressed  as: 


7+1 


D  2  2 

(7+3)  <£/or 

„  3/4 


+  4-7  (1  +  0.258  ~  )  . 


3/4 


<e> 

+  1.313  (  1/D)2  <7-l>  (1  +  3.67  =?  )  (-t— -  ) 

T0  <E),/8 


5.18 
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Substituting  the  relation  for  E  as  given  in  Equation  5.9  and 
assuming  a  Prandtl  number  of  0.72  the  final  form  of  the  drag  equation  is: 


^  +6.32^3,^, 


(1  +  0.258  rr  ) 
To 


M. 


(7  +3)  <  imy 

3/A 


(y+3)  [(  7-l)(l  +  2.6A  ^i)] 


1/8 


2  +  367  T5  > 

2.37  U/ny  (7  -1) 


M 


3/A  I  2 


[<7-l)(l  +  2.6A  \)\  1M  \  (R^ 


5.19 


which  can  also  be  written  in  the  form: 

H  ,  3/4 
C_ 


where , 


'M 


3 /A 


D 


(  1/D) 


2  +  a2 


V(V« 


+  a3  (  i  /D)‘ 


Wc 


a  ,40±  .1), 

1  (7+  3)2 


a2  '  6'37  I*  3) 


7+  1  'Ty 


(1  +  0.258  ^  ) 


(7-1) 


1/8 


(1  +  2.6A  jv  )1/8 

T  T0 

3/4  <l  +  3'67  T?  ) 

a  -  2.37  (7-l)J/4  - A-1/4 

J  (1  +  2.6A  J!)  ' 

T0 


5.20 


Equation  5.20  is  useful  in  making  quick  preliminary  drag  variation 
estimates  for  re-entry  trajectory  studies.  For  example,  assuming  a  7*  l.A 
and 


the  general  formula  for  all  cone  drags  (M»0C»1;  M>10;6C  <  2A°)  is: 


5-8 


AERONUTRONIC  DIVISION 


2 


„  .  0-496 

“D  (i/D)2 


Moo 


3/4 


Moo 


3/4 


+  2.88 


+  1.53  (£  /D) 


iVTvJ 


5.21 


Note  that  the  Inviscid  pressure  drag  term  above  (l.e.,  — ; r  )  Is  very 
nearly  that  given  by  Newtonian  theory  (l.e.,  (  £ /D) 


0.5 

(X/D)2 

The  principal  limitation  to  the  general  validity  of  Equation  5.20, 
aside  irom  the  approximations  discussed  in  the  text,  is  associated  with  the 
use  of  the  power  law  variation  of  viscosity  with  temperature. 

For  very  high  stagnation  temperatures  or  very  low  free-stream 
temperatures  (such  as  those  obtained  in  wind  tunnel  testing  at  high  Mach 
numbers)  the  power  law  variation  in  the  form  used  in  the  previous  deriva¬ 
tion  can  be  considerably  in  error.  The  simplified  form  of  the  drag  equation 
can  still  be  applied  with  considerable  success  if  a  suitable  transformation 
in  the  effective  Reynolds  number  is  made.  This  transformation  will  be  dis¬ 
cussed  further  in  Section  5.3  on  Correlation  of  Experimental  Data  and  con¬ 
sists  of  the  replacement  of  the  Reynolds  number  (RLoo)  in  Equation  5.20  by 
a  pseudo  or  effective  Reynolds  number  (R^  oo  eff  ective^  ^or  t*le  PurPoses 
the  drag  calculation.  To  anticipate  the  results  of  Section  5.3,  Equation 
5.33  has  been  utilized  to  determine  the  "effective"  Reynolds  number  and  a 
calculation  has  been  made  using  the  simplest  form  of  the  approximate  drag 
equation  (Equation  5.21)  applied  to  the  test  conditions  of  Reference  22. 

A  comparison  between  this  caluclation  and  the  experimental  data 
points  is  shown  in  Figure  5-1.  It  is  apparent  that  the  approximate 
theory  gives  very  good  results  at  this  high  Mach  number.  However,  since 
the  theory  is  based  on  hypersonic  approximations,  the  agreement  would  be 
less  favorable  at  lower  Mach  numbers. 

Finally,  It  must  be  remembered  that  interaction  theory,  whether 
exact  or  approximate,  is  no  longer  applicable  at  a  very  low  Reynolds 
number . 


5.2  DERIVATION  OF  CORRELATION  PARAMETERS  FOR  CONE  DRAG 

In  this  section  a  set  of  parameters  are  obtained  for  correlating 
the  drag  of  slender  cones  at  supersonic  and  hypersonic  speeds  . 
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FIGURE  5-1.  COMPARISON  BETWEEN  EXPERIMENT  AND  APPROXIMATE  THEORY 
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Ic  will  be  shown  that  through  the  use  of  these  parameters, 
experimental  data  for  a  wide  range  of  Mach  and  Reynolds  numbers  can  be 
unified  into  a  single  curve.  This  correlation  curve  shows  agreement  with 
interaction  theory  at  the  higher  Reynolds  numbers.  However,  at  very  low 
Reynolds  numbers,  significant  deviations  from  interaction  theory  are 
observed.  This  is  to  be  expected  since  it  is  well  known  that  as  the  flow 
density  decreases  continuum  flow  concepts  become  invalid  and  must  be 
replaced  by  those  appropriate  to  free  molecule  flow.  It  is  also  recog¬ 
nized  that  a  transition  regime  must  exist  between  the  continuum  and  free 
molecule  regimes. 

The  principal  advantage  of  the  present  correlation  procedure 
is  that  it  permits  the  consolidation  of  the  limited  available  low  density 
test  data  in  such  a  way  as  to  obtain  useful  trend  information  in  the 
transition  regime.  Such  information  is  particularly  useful  because,  at 
the  present  time,  there  is  no  well  established  theory  for  this  transition 
regime . 


The  correlation  parameters  to  be  discussed  here  are  based  on 
the  simplified  form  of  the  cone  drag  equation  (e.g.  Equation  5-20  of 
Section  5.1). 


It  is  observed  that  the  cone  drag  variation  Indicated  in  this 

Equation  is  characterized  by  three  different  regimes.  These  regimes 

are  related  to  the  numerical  magnitude  of  the  ratio  3/4. 

M°° 

v^r 

Li  00 


For  small  values  of  the  above  ratio,  the  second  and  third 
terms  of  the  drag  equation  are  negligible  compared  to  the  first  term, 
therefore,  the  cone  drag  is  simply  the  inviscid  pressure  drag.  For 
intermediate  values  of  the  ratio,  the  cone  drag  varies  as  a  function  of 
the  parameter  3/4  ,  while  for  large 


v<R rr 

Li  00 

values  of  the  ratio,  the 


drag  varies  as  a  function  of  M « 


3/2 


"l 


00 
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There  is  no  difficulty  in  obtaining  cone  inviscid  pressure 
drag  by  standard  methods.  Since  we  are  interested  principally  in  the 
rarified  flow  regimes  where  viscous  effects  are  important,  it  is  con- 
venient  to  eliminate  the  cone  pressure  drag  from  the  correlation  process. 
This  is  done  by  rewriting  Equation  5.20  in  the  form: 


+  a3  (  i/D)2 


5.22 


From  the  above  equation  it  is  obvious  that  we  could  plot  Cp-Cup 


3/4 


vs  a 


2  you 


and  obtain  a  general  correlation  for  small  to  moderate 


values  of  00  ,  however,  there  would  remain  a  divergence  at  larger 

values  of  this  ratio.  The  problem  is  to  extend  the  range  of  correlation. 
This  can  be  done  at  the  expense  of  a  more  complex  correlation  parameter . 


By  restricting  7  to  a  value  of  1.4  and  evaluating  the  coeffi¬ 
cients  &2  and  a3  f°r  various  values  of  Tw/T0,  it  is  found  that  &2  can  be 


T 

T 

approximated  as  82  *  2.82 

and  a^  «  1.12 

1+2.4  =* 

A0 

Substituting  these  values  into  Equation  5.22  and  neglecting  the 


temperature  contribution 

in  82  we  have: 

CD-°D 

P 

- 

2.62 

M  3/4 
n  00 

+  1.12 

1  +  2 

T 

.4  zr 

T 

0 

( <e/D)2 

1  mJ" 

5.23 
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This  suggests  a  correlation  parameter 


M 


3/4 


V  <*L> 


1  +  a  (i  /D)2  (1  +  b  ~ 

l0 


"oo 


3/4 


*V„ 


where  a  *>  0.4 

b. »  2.4 

This  Is  equivalent  to: 


where  c^,  the  pressure  coefficient,  is: 


5.24 


5.25 


c  .  t&JL  Jl_ 

P  ( 7+  3)2  (£/D)2 


2  sin2  9 

c 


The  above  representation  of  the  cone  pressure  Is  accurate  for 
very  high  Mach  numbers;  however,  at  lower  Mach  numbers  the  expression 
should  be: 


According  to  the  complete  Interaction  theory,  as  previously  given, 
P 

Q 

the  ratio  —  appears  directly  in  the  expression  for  the  cone  pressure 

“oo 

drag  (Section  2)  and  indirectly  In  the  expression  for  the  skin  friction 
(Section  4). 

In  the  correlation  procedure  being  used  here,  the  inviscid 
pressure  drag  is  subtracted  from  the  total  drag  at  the  outset;  therefore, 
assuming  this  inviscid  pressure  drag  is  correctly  evaluated  from  Kopal 
tables  or  similar  reference,  the  pressure  approximation  expressed  by 
Equation  5.25  does  not  influence  this  part  of  the  correlation. 
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Of  the  remaining  terms  associated  with  viscous  effects,  the 
largest  contribution  is  that  due  to  skin  friction.  From  Equation  4.6 
it  is  seen  that  the  principal  coefficient  of  the  skin  friction  equation 
is: 


cotQ  -  2  l  /D 


% 

‘loo 


553 


C 

c 


M 


M„ 


1/2 


Due  to  the  number  of  approximations  involved,  it  is  not  ex¬ 
pected  that  anything  more  than  the  form  of  the  correlation  parameter  X 
is  given  by  the  above  relation.  In  particular,  the  constants  a  and  b  are 
subject  to  adjustment . 

It  is  interesting  to  note  that  for  a  moderate  degree  of  flow 


rarification  the  appropriate  correlating  parameter  is  ■  . . . .  ,  rather 


than  the  usually  accepted  .  The  parameters  are  very  similar, 

YOU) 

Li  00 

however,  and  at  the  lower  Mach  numbers  the  difference  between  M  and 
M 

is  small;  hence  * ■=  as  will  provide  a  fair  correlation. 

V(V„ 

At  this  point  we  return  to  the  complete  interaction  theory  as 
given  in  Chapter  2-4.  If  the  parameter  X  is  to  correlate  experimental 
data,  it  should  also  correlate  points  calculated  by  the  complete  inter¬ 
action  theory.  This  is  certainly  true  for  the  intermediate  density  range 
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where  the  theory  is  known  to  agree  with  experimental  data.  It  is  also 
taken  to  be  a  requirement  in  the  low  density  regime  for  lack  of  a  better 
criterion. 


Figure  5-2  shows  the  correlation  obtained  for  theoretically 
calculated  drag  curves  assuming  various  combinations  of  free  stream  Mach 
number,  wall  temperature,  and  cone  angle.  The  correlation  parameter  X, 
as  shown  in  this  figure,  is  slightly  different  in  form  from  that  given 
in  Equation  5.24.  Aside  from  the  numerical  values  found  to  be  most 
suitable  for  the  constants  a  and  b  (a  ■  0.5,  b  *  2.5)  the  term 


i  +  LAM 

7  M  a, 


2 

2 


has  been  added  to  the  correlation  parameter. 


This  term  was 


found  necessary  to  obtain  agreement  at  the  lower  Mach  numbers  (M  <  15). 

It  is  basically  a  correction  to  the  skin  friction  contribution  to  account 
for  certain  pressure  effects  that  were  over-simplified  by  the  approxima¬ 
tions  used  in  the  derivation  of  Equation  5.20.  The  reasoning  behind  this 
correction  factor  is  discussed  in  the  following  text. 


recalled  that  in  the  derivation  of  Equation  5.20, 
number  was  assumed  hypersonic  so  that  the  invlscid 
could  be  represented  as: 

7  +  1  Moo2 

(  7  +  3  )  2  U/D)2 

we  can  write 

C 

P 

Q 

Using  the  hypersonic  approximation  for  —  ,  as  given  by  Equation  5.25 

*00 

we  have: 


It  will  be 
the  free  stream  Mach 
pressure  on  the  cone 


2  1/2  - - 

*  2  U D  (cp  )  'V Coo 


5.27 


5-15 


AEKONUTRONIC  DIVISION 


If,  on  the  other  hand,  we  uae  the  more  accurate  Equation  5.26  to  represent 
the  pressure  ratio,  we  have: 

_  ,  \l/2 


cot6 


hV 


Mo 


<*L> 


2  l  /D  (c 


1  + 


7  m  ( 


P  2 


1/2 


5.28 


It  is  apparent  that  Equation  5.28  is  identical  to  Equation  5.27 
except  for  t{ie  correction  factor 

,  1 1/2 


1  + 


C  'Y  H  00 
P 


Since  r  7  M  is  generally  small  compared  to  I,  the  correction 

factor  can  be  further  approximated: 

1/2 


1  + 


1 


7  M 


P  2 


1  + 


111. 


»  •  •  •  • 


P  2 


2  1 

or  expressing  as  2  sin  8  «  ^  ^ 


1  + 


1 


,  7  M°q 

'P  2 


1,2  «  ,  +  1  U.ffl)2  + 

7 


5.29 


This  is  the  form  of  the  correction  factor  that  must  be  applied 
to  the  correction  parameter  X  to  account  for  the  pressure  effects  on  skin 
friction. 
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It  was  found  in  the  process  of  correlating  the  theoretically 
calculated  drag  curves,  that  best  agreement  was  obtained  when  the  constant 
2  was  dropped  from  the  correction  factor.  This  is  probably  due  to  addi¬ 
tional  implicit  approximations  which  have  not  been  directly  accounted  for. 
(The  assumption  of  a  constant  value  for  the  coefficient  F;  for  example.) 

In  final  form  the  correlation  parameter  is: 


M  oo 


3/4 


1  +  (  l/D Y 
2 


1  +  2.5  Jlw 

To 


M® 


3/4 


1  +  LLL pi 

7  Moo2 


5.30 


It  is  apparent  from  Figure  5-2  that  good  correlation  i3  obtained 
between  the  various  theoretically  calculated  drag  curves. 

We  now  proceed  to  the  correlation  of  experimental  data. 


5.3  CORRELATION  OF  EXPERIMENTAL  DATA 

Most  of  the  available  hypersonic  experimental  test  data  have 
been  obtained  at  very  low  freestream  temperatures  (i.e.,  50°  to  100°K). 
This  requires  that  some  caution  be  exercised  in  use  of  the  correlation 
parameter  X  as  defined  in  the  previous  section. 

The  problem  is  centered  about  the  fact  that  the  viscous  drag  of 
a  cone  is  a  function  of  local  conditions  behind  the  shock  while  the  corre¬ 
lation  parameter  X  is  based  on  conditions  in  front  of  the  shock.  In  the 
derivation  of  X,  the  approximate  equation. 


<Rx>  P  I  M  \  1/2 

<»„/  '  p^T  ( ]  5- 

was  used  to  express  the  Reynolds  number  transformation  across  the  shock. 
This  equation  is  based  on  the  power  law  temperature-viscosity  relationship 
(  O)  ■  3/4)  and  assumes  that  the  Reynolds  number  ratio  is  Independent  of  the 
freestream  temperature.  This  is  not  strictly  true,  particularly  at  the 
lower  freestream  temperatures.  A  more  accurate  representation  of  the  Rey¬ 
nolds  number  ratio  across  the  shock  is  given  by  the  equation  (based  on  the 
Southerland  relation): 
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(R  ) 

X 

c 

p 

c 

M 

c 

T  00 

ip 

[  1  ♦  “°°K 

1 00 

Too 

T 

c 

<V  ' 

A  00 

Poo 

Moo 

l 

c 

1  +  110°K 

T 

L.  *■  oo 

* 

According  to  this  latter  representation,  a  given  viscous  drag 
level  (which  corresponds  to  a  fixed  unit  Reynolds  number  behind  the 
shock)  occurs  at  different  values  of  the  freestream  unit  Reynolds  number 
as  the  freest ream  temperature  Is  changed. 

Alternately,  for  a  fixed  freestream  Reynolds  number,  the  cone 
drag  varies  as  a  function  of  freestream  temperature. 

Because  the  correlation  parameter  X  Is  based  on  the  approximate 
Equation  S.31,  this  variation  with  freestream  temperature  Is  not  accounted 
for. 


It  Is  noted,  however,  that  if  an  effective  freestream  Reynolds 
number  Is  substituted,  so  that  the  power  law  gives  the  correct  downstream 
conditions,  then  the  correlation  parameter  X  Is  valid. 

The  relationship  between  the  effective  Reynolds  number  that 
must  be  used  for  X  and  the  true  freestream  Reynolds  number  as  measured 
In  a  test  is: 


!1+  noV 

Too 

1  T 

T 

m 

Effective  Measured 

l  OO  1 

c 

'.  7  1 1 0°K  1 

1  +  T 

Aoo 

Uc  1 

5.33 


It  may  be  observed  that  as  the  freestream  temperature  approaches 
0  the  above  equation  reduces  to 


0; 


effective 

*  j 

measured 


3/4 
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while  for  large  freestream  temperatures: 


T 


00 


^effective 

^measured 


1/4 


Between  these  two  limits  there  must  be  a  temperature  range  where 
R 


R 


effective 

mensured 


1.  This  is  the  temperature  range  where  the  power  law 


viscosity  relationship  is  accurate.  For  slender  bodies  this  usually  occurs 
between  200°  and  300°K. 


The  corrected  correlation  parameter  which  units  the  experimental 
data  is,  then: 


The  last  obstacle  to  be  overcome  in  order  to  achieve  correlation 
of  the  experimental  data  is:  The  inviscid  cone  pressure  drag  level  Cpp, 

is  different  for  different  cone  angles  and  Mach  numbers.  Therefore,  even 
assuming  that  the  cones  approach  a  uniform  drag  level  (i.e.,  -f  * 

2.0),  as  the  Reynolds  number  goes  to  zero,  the  asymptote  of  the  variable 

Cn  -  Cn  will  lie  at  different  levels,  and  correlation  at  the  lower 

Dexp  DP 

Reynolds  numbers  will  not  be  obtained.  This  difficulty  can  be  avoided  by 
using  the  quantity  (Cp  -  Cp  )/(Cp  -  Cp  )  as  the  ordinate  parameter.  At 

the  same  time,  a  modification  must  be  made  to  the  abscissa  to  retain  corre- 
lation  at  high  Reynolds  numbers.  This  modification  consists  of  replacing 
X  by  X/Xpjj  ,  where  Xpjj  is  the  value  of  X  at  the  point  where  the  interaction 
theory  drag  curve  passes  through  the  free  molecule  drag  level. 
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In  practice  this  reduces  to: 

CD  FM  ‘  CD_ 

XFM-  - 2T5 -  5-35 


Figure  5.3  shows  the  correlation  between  the  experimental  drag 
data  obtained  for  various  cones  at  various  Mach  numbers,  Reynolds  numbers, 
and  wall  temperature  ratios.  It  is  apparent  that  excellent  agreement  is 
obtained  between  the  various  test  conditions  and  that  a  single  curve 
representing  the  viscous  interaction  drag  can  be  faired  through  all  the 
data  points.  The  particular  shape  of  the  faired  curve  represents  addi¬ 
tional  information  not  presented  here.  This  information  is  related  to 
sphere  drag  correlations  which  will  be  the  subject  of  a  separate  memorandum. 

With  the  correlation  curve  presented  here,  the  drag  of  a  given 
cone  under  given  conditions  can  be  calculated  by  evaluating  X/Xp^  for  the 

appropriate  values  of  M,  R  and  TwTq  and  reading  the  corresponding  value  of 


The  drag  calculated  in  this  manner  is  not  restricted  to  Rey¬ 
nolds  number  ranges  where  interaction  theory  is  known  to  apply.  Devia¬ 
tions  from  interaction  theory  due  to  slip  flow,  etc.  are  accurately  re¬ 
flected  and  extrapolations  into  the  transition  regime  can  be  made  with 
confidence . 
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CHAPTER  6 

FLOW  OVER  CONES  AT  ANGLE  OF  ATTACK 


6.1  GENERAL 

The  complexity  of  the  three-dimensional  flow  fields  occurring 
over  cones  at  arbitrary  angle  of  attack  is  so  great  that  rigorous  solu¬ 
tions  to  the  equations  of  motion  have  not  yet  been  obtained.  This  is 
not  a  deterrent  to  many  engineering  designs  since  certain  restrictive 
assumptions  can  be  made  and  solutions  obtained  for  special  cases.  For 
example,  in  hypersonic  flow,  the  assumption  of  high  Reynolds  numbers 
permits  an  estimate  of  the  aerodynamic  forces  on  bodies  by  neglecting 
skin  friction  effects.  The  necessity  of  including  the  skin  friction 
forces  in  low-density  flows  due  to  boundary  layer  effects,  however,  is 
an  added  complication  which  requires  considerably  more  effort  to  solve. 
When  these  boundary  layer  effects  become  important,  the  determination 
of  the  flow  field  must  consider  not  only  cross  flow  and  real  gas  effects, 
for  sufficiently  high  angle  of  attack,  but  also  the  viscous  interaction 
occurring  at  the  smaller  angles  of  attack.  This  latter  phenomenon 
causes  a  local  increase  in  the  magnitude  of  both  pressure  and  skin 
friction  forces  as  well  as  displacing  the  shock  wave  to  perturb  the 
local  temperature-density  field. 

If  the  freestream  density  is  sufficiently  low,  the  inclusion  of 
the  skin  friction  forces  may  still  be  required  at  high  angles  of  attack, 
even  though  any  viscous  interaction  has  disappeared  due  to  an  effective 
blunt  body  flow  situation.  The  determination  of  the  flow  field  at  inter¬ 
mediate  angles  of  attack  may  be  subject  to  serious  errors  since  the  local 
flow  properties  must  be  determined  accurately  to  estimate  the  magnitude 
of  the  low-density  flow  effects. 
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The  question  then  naturally  arises  as  to  what  would  be  the  proper 
approach  in  examining  the  flow  field  over  slender  cones  at  arbitrary  angle 

of  attack  for  current  engineering  designs  is  controlled  by  the  following: 

If  primary  consideration  is  given  to  the  aerodynamic  forces  acting  on  the 
body,  then  a  detailed  examination  of  the  distribution  of  local  flow  vari¬ 
ables  may  not  be  required.  In  the  present  case,  primary  emphasis  will  be 
placed  on  the  estimation  of  the  aerodynamic  force  coefficients.  An  attempt 
will  be  made  to  estimate  local  flow  properties  from  a  brief  examination  of 
the  shock  shape  insofar  as  these  are  used  to  predict  the  force  coefficients. 
The  general  approach  will  be  to  examine  two  limiting  cases:  (1)  the  small 
angle  of  attack  case,  where  cross  flow  effects  are  small  or  negligible,  and 
(2)  the  high  angle  of  attack  case,  where  cross  flow  predominates  and  any 
viscous  interaction  effects  are  negligible. 

The  small  angle  of  attack  case  may  be  considered  to  be  a  pertur¬ 
bation  of  the  results  already  presented.  It  will  be  shown  that  the  appli¬ 

cability  of  the  predictions  extend  over  a  range  of  angle  of  attack  up  to 
the  order  of  the  cone  semi-vertex  angle.  Progressive  increases  in  the 
angle  of  attack  beyond  this  point  produces  increasing  error  as  the  cross 
flow  increases.  The  cone  at  high  angle  of  attack  (90  >a»0)  will  be  ex¬ 
amined  by  the  application  of  a  modified  infinite  swept  cylinder  theory. 

This  greatly  simplifies  the  problem  since  similar  solutions  are  available 
to  the  boundary  layer  equations  for  swept  cylinders,  which  include  variable 
properties  and  the  inter-dependence  of  the  axial-and-cross-f low.  This 
approach  has  previously  been  used  in  calculating  the  stagnation  line  heat 
transfer  on  a  slender  cone  .  In  the  intermediate  angle  of  attack  range, 
the  assumption  will  be  made  that  a  smooth  interpolation  curve  between  the 
limiting  cases  is  adequate.  This  simple  patching  procedure  will  be  made 
by  arbitrarily  fairing  the  curve  tangent  to  both  calculated  portions. 


6.2  LAMINAR  BOUNDARY  LAYER  ON  A  CONE  AT  ANGLE  OF  ATTACK. 

The  laminar  boundary  layer  flow  about  a  circular  cone  at  small 
angle  of  attack  has  been  analyzed  in  the  plane  of  symmetry  by  several  in¬ 
vestigators^  *30,3^ .  Each  points  out  the  difficulty  of  obtaining  solutions 
away  from  the  stagnation  line.  Moore30,  however,  has  a  discussion  on  the 
general  behavior  of  the  boundary  layer  over  a  cone  at  angle  of  attack 
which  is  pertinent  to  the  present  problem.  He  shows  that  the  boundary 
layer  tends  to  be  thicker  on  the  top  of  the  cone  at  angle  of  attack. 

Quoting  from  Moore,  "This  is  to  be  expected  since  the  fluid  near  the  base 
of  the  boundary  layer  has  low  inertia,  is  therefore  inclined  to  follow 
the  direction  of  the  circumferential  pressure  gradient  more  closely  than 
in  the  outer  flow,  and  thus  tends  to  drain  away  from  beneath  the  cone  and 
accumulate  near  the  top.  No  separation  is  encountered  because,  for  small 
angle  of  attack,  the  pressure  gradient  is  always  favorable*". 


•ff 

Approximately  for  a  <  cone  half-angle,  0. 
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"For  larger  angles  of  attack,  when  the  pressure  gradient  reverses 
direction  near  the  top  of  the  cone,  experiment  indicates  the  formation  of 
boundary- layer  'lobes'.  When  the  angle  of  attack  is  further  increased  the 
lobe  pattern  finally  breaks  away  from  the  body  to  form  a  vortex  street." 
This  may  be  illustrated  by  the  following  cross-sectional  views  of  the  boun¬ 
dary  layer  on  a  cone  at  various  angles  of  attack: 


The  angle  of  attack  where  these  different  conditions  occur  may 
be  considerably  modified  in  the  present  case  by  the  very  high  Mach  numbers 
and  effects  of  viscous  interaction  on  slender  cones.  In  the  limiting  case 
of  small  angle  of  attack,  a  major  consideration  of  the  analysis  is  the  de¬ 
termination  of  what  angle  of  attack  does  crossflow  become  non-negligible 
in  estimating  the  aerodynamic  forces.  If  streamlines  are  traced  over  a 
conical  body  at  angle  of  attack  by  following  the  pressure  gradient,  the 
turning  angle  of  the  flow  (which  is  a  measurement  of  the  crossflow)  in¬ 
creases  with  increasing  distance  away  from  the  meridian  plane  of  symmetry 
or  stagnation  line.  A  rigorous  analysis  to  determine  the  magnitude  of  the 
crossflow  should  therefore  include  the  entire  surface  of  the  cone.  This 
cannot  be  accomplished  for  arbitrary  a.  A  simple  estimate  of  the  cross- 
flow  magnitude,  however,  may  be  made  for  the  present  analysis  in  the  fol¬ 
lowing  manner: 

(1)  A  criteria  for  negligible  crossflow  is  arbitrarily 
assumed  by  taking  the  crossflow  velocity  gradient 
at  the  stagnation  line  as  zero.  Solutions  are  avail¬ 
able  for  this  case;  at  least  in  the  case  where  the 
boundary  layer  is  thin  and  viscous  interaction 
unimportant . 


6-3 


AERONUTRONIC  DIVISION 


(2)  The  analysis  will  estimate  the  magnitude  of  the 
stagnation  line  crossflow  velocity  gradient  as 
a  function  of  freestream  Mach  number,  cone  half- 
angle  and  a . 

We  have  already  confined  ourselves  to  the  case  of  M  >>  1  and 

oo  00 

slender  cones  (i.e.,  20  >a>5  ).  The  validity  of  the  estimate  for  neg¬ 
ligible  crossflow  will  be  subsequently  determined  by  comparison  with  ex¬ 
periment  insofar  as  estimates  of  aerodynamic  force  coefficients  are 
concerned.  The  crossflow  velocity  gradient,  W^,  at  the  stagnation  line 
has  been  calculated  in  the  MIT  tables^  ^  and  reproduced  in  Reference  15 
for  a  range  of  conditions.  These  calculations  Indicate  a  linear  decrease 
in  with  increasing  Mach  number  at  a  given  value  of  a  and  0.  The  ex¬ 

amination  of  5-  and  10-degree  half-angle  cone  values  of  presented  in 
Reference  15  further  indicates  that  would  be  very  small  or  negligible 
for  a  of  the  order  of  or  less  than  a  .  This  will  then  be  taken  as  the 
limiting  region  of  small  angle  of  attack.  This  region  will  be  treated 
by  assuming  negligible  crossflow,  i.e.,  assuming  the  boundary  layer  grows 
from  the  cone  vertex  as  in  the  case  a  *  0  providing  local  properties  are 
used  at  each  point  in  the  flow  field.  The  results  of  the  zero  angle  of 
attack  case  may  be  utilized,  assuming  that  the  local  shock  slopes  are  in¬ 
dependently  known.  It  is  to  be  expected  that  the  increasing  error  occurr¬ 
ing  with  increasing  angle  of  attack  is  perhaps  not  as  large  or  important 
in  the  determination  of  aerodynamic  force  coefficients  where  local  property 
variations  are  integrated  out,  as  in  the  determination  of  local  heat  trans¬ 
fer  coefficients  which  depend  directly  on  local  phenomena. 

As  previously  noted,  solutions  for  the  other  limiting  region  of 
large  angle  of  attack  will  be  developed  from  infinite  swept  cylinder  theory 
for  the  slender  cones  and  neglecting  the  effects  of  any  viscous  interaction. 
Beginning  at  a  *  90  degrees,  the  solution  is  a  two-dimensional  blunt  body 
problem.  As  a  i8  decreased,  the  contribution  of  skin  friction  to  the  forces 
becomes  more  important,  at  least  in  the  axial  direction.  The  question  as 
to  the  region  of  importance  of  viscous  interaction  then  arises.  Calculations 
have  been  made  of  the  interaction  parameter,  *  0,  for  0*5  and  10 
degrees  and  as  a  function  of  altitude  in  Figure  6-1  for  »  20  and  x  *  1  ft. 
The  method  for  obtaining  local  shock  slope  at  the  stagnation  line  (and  re¬ 
sultant  flow  properties)  will  be  examined  in  Section  6.3.  Figure  6-1  indi¬ 
cates  that  for  the  slender  cones  of  interest  fairly  definitive  limits  may 
be  set  for  the  lowest  angle  of  attack  at  which  viscous  interaction  effects 
should  not  appreciably  influence  the  forces.  For  example,  based  on  the 
interaction  regions  already  established  (criteria  here  being  taken  as  limit¬ 
ing  value  of  y  ■  one-tenth  of  the  value  of  X  for  onset  of  strong  interaction), 
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an  insulated  10-degree  half-angle  cone  should  ^rtot  have  any  appreciable 
viscous  interaction  effects  noticeable  for  a  >  45  degrees  at  200  kilo- 
feet  altitude;  whereas  this  would  change  for  a  highly  cooled  cone  to 
a  10  degrees  for  the  same  conditions.  The  stagnation  line  region  just 
examined  represents  the  region  of  maximum  compression  at  angle  of  attack, 
and  consequently,  the  region  of  lowest  viscous  interaction  on  the  surface 
of  the  cone.  This  is  particularly  true  at  sufficiently  low  angles  of 
attack  where  the  flow  originates  at  the  cone  vertex,  instead  of  at  the 
stagnation  line. 


6.3  SHOCK  SHAPE  AT  ARBITRARY  ANGLE  OF  ATTACK 

In  Section  6.2,  requirements  were  outlined  for  obtaining  local 
flow  properties.  A  considerable  knowledge  of  the  properties  in  the  hyper¬ 
sonic  flow  field  may  be  obtained  from  the  local  slope  of  the  shock  wave. 

The  conical  shock  angle  at  a  ■  0  was  determined  numerically  in  Reference  32 
and  graphically  presented  in  NACA  Number  TR1135  for  a  perfect  gas.  The  shock 
shape  was  also  obtained  in  References  33  and  34  for  very  small  angles  of 
attack;  but  these  are  insufficient  for  the  present  problem.  A  past  for  an 
approximate  estimation  of  the  shock  angle  in  the  plane  of  symmetry  is  to 
use  an  equivalent  cone  method  where  the  shock  angle  is  obtained  from  an 
equivalent  cone  with  a  equiv  *  a  +  a  .  This  approximation  becomes  pro¬ 
gressively  worse  with  increasing  a  ,  and  approaches  the  two-dimensional 
wedge  value  at  large  a-  The  use  of  hypersonic  small  perturbation  theory 
is  also  inadequate,  as  the  assumptions  for  its  validity  are  violated  ex¬ 
cept  at  very  small  angle  of  attack. 

35  36 

More  recently,  Gonor  ’  has  examined  the  conical  flow  field 
at  angle  of  attack  for  M  >>1,  and  derived  an  analytical  expression  for 
the  shock  angle  in  the  plane  of  symmetry.  This  expression  has  the  form: 

1  + - 2~  2 -  I  6A 

(  7  -ljM^  sin  (a  +  a)J 

where  t  »  +  sin  a/sin  a  cos  (  a  ±  a  ) 


a  + 


tan(  pig)  [ 

t2  L 


(l+t)  In  (1+t) 
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The  upper  sign  In  Equation  6.1  refers  to  a  meridional  plane  located  on  the 
most  windward  side,  and  the  lower  sign  to  a  meridional  plane  on  the  most 

2 


leeward  side.  The  term 


7- 

7+1 

across 


( 7  -1)M2  sin2(a  +  a)  _ 
the  conical  shock  wave  in 


is  approximately 
a  very  high  Mach 


the  density  ratio,  e 
number  flow.  Gonor  has  compared  Equation  6.1  with  experiment  and  found 
satisfactory  agreement  at  small  angles  of  attack  at  Mffl  ■  4.0.  The  leeward 


side  solution  is  invalid  for  a  ><j*  the  windward  side  solution  has  been 
examined  for  a>Cand  is  compared  with  the  equivalent  cone  method  in 
Figure  6-2  for  c  ■  10  degrees  at  M^*  ».  Both  solutions  exhibit  sudden 
asymtotlc  increases  as  a  increases.  In  the  region  of  large  a,  9g  should 

approach  an  equivalent  value  for  a  blunt  body  standoff  distance.  The  limit¬ 
ing  value  of  6a  at  large  a,  for  which  Equation  6.1  can  be  expected  to  be 
valid,  is  uncertain.  Such  experiments  as  those  in  Reference  37,  however, 
indicate  a  near  conical  flow  up  to  an  angle  of  attack  of  approximately 
60  degrees . 


An  estimation  of  the  shock  angle’s  limiting  value  at  very  large 
angle  of  attack  may  be  made  by  considering  the  flow  over  a  swept  cylinder 
at  zero  yaw.  The  shock  standoff  distance  varies  inversely  as  the  local 
radius  of  curvature.  If  this  type  solution  is  applied  to  a  slender  cone 
by  using  an  equivalent  local  radius,  a  shock  standoff  distance  is  gener¬ 
ated  which  varies  linearly  with  axial  distance.  It  may  be  considered  to 
be  a  shock  angle  in  the  conical  sense.  A  solution  of  this  type  will  now 


In  the  above  sketch,  the  cone  is  assumed  to  act  as  a  cylinder 
whose  local  radius  of  curvature,  R,  ,  is  parallel  to  the  flow  normal  to 
the  shock  wave.  D 
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Then 

r.  “x  tan  0 

D 

“  rb  *l.n(90-  0)/sin(  a  +  a  +  A  ) 
-  x  sin  0 / cos(  0  -  0  ) 
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l 

sin  6  cos(  0  -  0)  -  (1+  — -  )sin  0 
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6.3a 

6.4 
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6.6 


Equation  6.6  is  an  implicit  relationship  for  0 g  which  must  be  solved  by 
trial  and  error  or  an  iteration  procedure.  It  is  seen  that  in  a  hyper¬ 
sonic  flow,  cos  (  0  -  0 )  %  1,  so  that  a  good  approximation  for  Equation 

6.6  is:  8 


A 

sin  6  ~  (1  +  - — )  sin  0  6  •  6a 

8  Rb 

This  solution  for  the  shock  angle  at  very  large  a  is  very  similar  to  that 
developed  in  Reference  38.  An  empirical  relationship  for  the  shock  stand¬ 
off  distance,  A,  is  also  presented  in  Reference  38: 


2  , \2 


“  0.38  +  2.4/ (Mjj  -1)  +  1.4/ (Mj,  -1) 

Substitution  of  6.7  into  6.6  gives 

sin  0scos(  0s-  0)  -  [l. 38+2.4/ (M^-l)  +  ^/(M^-l)2  ]  sin  0 
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or 


sin  0  % 


2.4  1.4 

1.38  +  — * -  + 


2  2  2 
**  "  1  (mn  - 


sin  0 
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The  use  of  equation  6.8  in  connection  with  the  method  of  Gonor 
(Equation  6.3)  has  been  compared  with  experimental  values^®  at  3.86 

in  Figure  6-3.  The  agreement  generally  good,  except  at  a  -  0  for  the  5- 
degree  cone.  Calculations  show  that  values  obtained  by  Gonor's  method 
are  in  poor  agreement  with  NACA  Number  1135  at  a  -  0  for  tan  a  1 

At  a  -  0,  values  from  NACA  Number  1135  should  therefore  be  used.  At 
large  angles  of  attack,  the  procedure  of  fairing  a  smooth  curve  tangent 
to  the  Gonor  curve  and  through  the  point  at  90°  -  Qg  obtained  from  the 
Mjj  method  produces  good  agreement  with  the  experimental  data.  Results 

for  a  -  5,  10  and  20  degrees  are  presented  in  Figure  6-4  for  M  -  10, 

25  and  infinity.  00 


At  small  a  the  shock  wave  will  still  be  approximately  conical 
in  shape,  although  the  shock  angle  on  the  windward  side  may  be  greater 
than  on  the  leeward  side^.  A  cross-section  of  the  shock  wave  tends  to 
be  elliptical  in  shape.  An  estimation  of  the  local  shock  angle  at  a 
point  exactly  90  degrees  away  from  the  plane  of  symmetry  gives  an  approx 
imate  value  of  that  at  a  *  0,  based  on  the  local  body  inclination  to 
the  freestream.  The  local  inclination  of  a  6  -  90  degree  generator  is 


where 


effective 


6. 


Calculations  for  a  ■  5  and  10  degrees  show  that  o  effective  to  the 
0-90  generator  decreases  very  slightly  at  small  a  : 
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FIGURE  6-4.  SHOCK  ANGLE  ON  CONE  AT  ANGLE  OF  ATTACK 
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Within  the  accuracy  of  .he  negligible  cross  flow  approximation  and  the 
distribution  assumptions  to  be  subsequently  employed  in  estimating  the 
aerodynamic  force  coefficients,  it  is  felt  that  0-90  degrees  corres¬ 
ponds  to  a  “0  for  a  reasonable  engineering  approximation. 

At  large  a ,  the  shock  angle  in  the  plane  of  symmetry  has 
been  discussed.  The  need  to  know  the  shock  slope  in  order  to  determine 
local  properties  is  not  so  great  as  in  the  case  of  small  a.  At  large 
a  ,  a  major  portion  of  the  flow  has  passed  through  the  shock  wave  near 
the  plane  of  symmetry,  and  then  mixes  with  the  flow  coming  through  the 
local  portion  of  the  shock  wave  as  the  flow  expands  around  the  body.  An 
expanding  blunt  body  flow  has  been  previously  examined^  and  the  results 
are  available  from  the  author.  The  detached  shock  shape  for  cylinders 
has  been  estimated  by  Love  .  The  pertinent  equations  as  applied  to 
the  case  of  a  slender  cone  at  large  angle  of  attack  are  presented  in 
Appendix  A. 
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CHAPTER  7 

THEORETICAL  ESTIMATES  OF  AERODYNAMIC  FORCE  COEFFICIENTS 
AT  ANGLE  OF  ATTACK 


7 . 1  AXIAL  FORCES  AT  SMALL  ANGLE  OF  ATTACK 

The  major  assumption  to  be  employed  in  analyses  of  the  present 
kind  is  the  one  with  respect  to  what  are  the  distributions  of  pressures 
and  skin  friction  shear  forces  over  the  surface  of  the  cone.  This  con¬ 
sideration  normally  leads  one  to  the  conclusion  that  the  solution  may  be 
divided  into  regions,  a  <  a,  and  a  >  u  (where  the  back  side  of  the  cone 
is  in  shadow).  In  the  second  region,  the  back  side  pressures  are  very  low, 
and  the  usual  procedure  is  to  neglect  it  in  any  consideration  of  forces 
acting  on  the  body.  A  further  division  of  the  solution  can  be  made  by 
examining  the  front  half  (-it/2  <  ♦  <  jt/2 )  and  back  half  (x/2  <  ♦  <  3rt /2 ) 
of  the  cone  separately.  In  low  density  flows,  however,  the  effects  of 
viscous  interaction  complicate  this  picture.  As  pointed  out  in  Chapter  2, 
the  magnitude  of  these  effects  increase  as  one  proceeds  away  from  the 
compression  region  at  angle  of  attack  towards  the  back  side.  It  is  not 
immediately  apparent  what  the  combined  effects  of  low  pressures  (or  local 
dynamic  pressure)  and  large  viscous  interaction  effects  are  on  the  back 
side,  and  at  what  angle  of  attack  the  back  side  can  finally  be  ignored. 

It  is  logically  assumed  that  the  pressure  contribution  for  a  >  a  to  the 
axial  force  would  be  small  although  not  necessarily  negligible  since  the 
trigonometric  factor  for  determining  the  pressure  component  varies  as 
sin  a.  The  same  thing  cannot  be  said  for  the  skin  friction  contribution, 
where  the  trigonometric  factor  is  cos  a-  This  problem  will  be  examined 
further  in  Chapter  7.6  when  comparisons  of  theory  and  experiment  are  made. 

In  the  determination  of  the  axial  force  coefficient  the  follow¬ 
ing  assumptions  will  be  made: 
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(1)  Modified  Newtonian  pressure  distribution  for  the 
lnviscid  pressures, 


-£—  =  cos2  0 


p*=l 


4>=0 


cos  0  =  (sin  cr  cos  0!  +  cos  a  sin  a  cos  ♦)/sin  (a  +  a) 


♦=0 


to  be  determined  from  the  shock  slope. 


(2)  Induced  pressures, 


p - 1  +  Fc  Xc  (Reference  Equation  2.22.) 


(3)  Linear  variation  of  viscous  interaction  effects  for 

0  <_♦  <  it/2,  and  it/ 2  <  ♦  <  n  as  determined  by  the  inter¬ 
action  parameter,  xc .  Also,  a  linear  variation  of  local  dynamic  pressure 
and  skin  friction  coefficient  with  $. 

The  pressure  contribution  to  the  axial  force  coefficient  is 
assumed  to  be 


A  -  q  A  J  PAdA  =  q  A  .  [  P  si"2  ad*xdx 

-  ref  x  ref  -*0  0 


cp  =  1 


With  assumptions _1  to  3,  Equation  7.1  becomes 


cf  -  2 
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sin  od  $  xdx 
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The  first  integral  represents  the  contribution  of  the  front  side  of  the 
cone,  and  the  second  integral  that  of  the  back  side.  Integration  of  7.2 
results  in,  for  the  front  side, 


For  the  back  side,  the  integral  will  be  assumed  to  hold  for  a  <  ct: 


Similarly,  the  skin  friction  contribution  to  the  axial  forces  is 


CA  =  o — A -  I  Cf(x>  cos  a  dA  7.5 

»  ref  -'A 

The  value  for  Cf(x)  and  resultant  integration  was  taken  from  Equations  4.2 
and  4.3,  where  both  the  viscous  interaction  and  transverse  curvature  effects 
are  included: 
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The  subscripts  0,  «/2,  n  refer  to  ♦  ,  and  the  values  of  (Cp*/2)  q^/qg  may  be 
obtained  from  Equation  4.3  by  using  local  properties  at  each  value  of  ♦  . 

The  total  axial  force  coefficient  is 


C 


T 

A 
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where  any  base  contribution  has  been  neglected.  Any  base  pressure  would 
subtract  from  the  value  of  given  in  Equation  7.9. 


7.2  NORMAL  AND  MOMENT  FORCES  AT  SMALL  ANGLE  OF  ATTACK 

The  same  assumptions  will  be  employed  for  the  normal  and  moment 
forces  as  in  the  determination  of  the  axial  forces.  The  integral  for  the 
normal  force  coefficient  then  becomes 

2  rL  c*  f 

C  =  - - -  /  I  p  cos  g  Sin  g  cos  odtxdx  - 

N  %  Aref  Jo  Jo  1  N 

Cj(x,  *)  q^($)  cos  ♦  sin2  adGxdxl  7.9 


Integration  of  Equation  7.9  gives 
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Equations  7.11  and  7.13  are  assumed  valid  for  a  <  a.  By  apply¬ 
ing  the  same  assumptions,  the  Integral  for  the  moment  coefficient  takes 
the  form: 
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A  positive  moment  is  taken  here  as  the  clockwise  direction. 
The  integration  of  Equation  7.15  gives 
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where 


Hx  -  (4*  -  8)  q#  +  (8  -  2k)  q^ 
♦  2 


H2  -  (8  -  2k)  q#  +  (K  -  8)  qn 
*  2 


-(1+  FcX^)3/2a^- <0.266  +  0.332?  0.517  +  0.913  — + 


0.121(y  -  1)  M‘j  -  0.232  +  FX* 

1.660  ?  jo. 517  +  0.913  ^  +  0.121  (y  -  1)  M*j 


1.162  F^ 


FCX*  YFCX, 

-f"  +  (1  +  W 


—  -1/2  1  I 

+  F  X J  '  +  yj-  - 


(1+  FcX,>1/2, 


in 


i  +  a  +  Fcy1/Z 
<  Fcx*>l/2 


=  sin  ct  cos  a/ sin  (a  +  Cl) 
=  cos  a  sin  a/sin  (a  +  a) 


T  P 

r*  *  n* 


°M  ‘  CM1  +  +  \  +  \ 


7.20 


7.3  SMALL  ANGLE  OF  ATTACK,  a  >  a 

In  Equations  7.4  to  7.11  the  solutions  for  the  back  side  of  the 
cone  apply  only  for  a  <  a.  In  the  region  a  >  a,  the  determination  of  the 
terms  containing  yFcX*  poses  a  problem  as  the  original  derivations  are 
valid  for  >  1.  Since  the  contribution  of  the  back  side  is 
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expected  to  be  smaller  than  that  of  the  front  side,  the  determination  of 
the  viscous  Interaction  effects  will  be  made  by  assuming  a  linear  variation 
of  Fc/y  =  0.492  at  Kb  =  1  to  Fc/y  =  1  at  Kfa  =  0,  in  a  manner  analogous  to 
the  case  of  the  flat  plate.  In  the  region  a  >  cr,  but  still  at  small  a, 
the  determination  of  the  back  side  contribution  to  the  force  coefficients 
is  liable  to  much  uncertainty.  The  general  level  of  the  back  side  pressure 
is  probably  of  the  order  of  the  free  stream  value.  The  consideration  of 
the  back  side  contribution  would  be  approximately  equivalent  to  the  use 
of  a  pressure  differential  for  the  pressure  level  on  the  front  side  in 
the  force  integrals  (Ca/A(P  -  Pw)dA)  and  neglecting  the  back  side  alto¬ 
gether,  if  there  were  no  viscous  interaction  effects. 


It  is  obvious  that  the  general  trend  of  the  back  side  contribu¬ 
tion  is  to  decrease  with  increasing  angle  of  attack.  The  pressure 
coefficient  across  a  shock  wave  varies  as  the  square  of  the  sine  of  the 
local  shock  angle.  It  will,  therefore,  be  assumed  that  the  back  side 
pressure  is  reduced  approximately  by  cos2  a.  An  estimation  of  the  back 
side  contribution  of  skin  friction  forces  will  then  be  made  by  assuming 
freestream  conditions  at  ♦  =  jt.  a  =  a,  and  allowing  the  skin  friction 


coefficient  to  decrease  as  cos?  at  so  that 
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cos2  a  in  the  case 


a=0 

of  the  axial  skin  friction  force  coefficient.  An  examination  of  the  skin 
friction  forces  contributing  to  and  Cy  reveals  the  front  and  back  side 
are  of  the  same  order  of  magnitude  and  of  opposite  sign.  The  net  contri¬ 
bution  of  shearing  forces  is  negligible.  A  similar  estimation  in  the 
case  of  the  back  side  pressure  contribution  will  be  made  from  the  a  =  a 
condition  by  assuming  a  constant  inviscid  pressure  on  the  back  side  of 
freestream  magnitude.  This  gives  for  a  >  cr : 
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where  the  interaction  parameters  are  to  be  evaluated  at  a  =  a. 
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7.4  CENTER  OF  PRESSURE 


The  center  of  pressure 
previous  normal  and  moment  force 


for  the  cone  may  be  determined  once  the 
coefficients  are  known.  It  is  given  by 
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7.5  AERODYNAMIC  FORCES  AT  LARGE  ANGLE  OF  ATTACK 


In  the  large  angle  of  attack  region,  the  viscous  interaction 
effects  have  essentially  disappeared.  In  addition,  any  contribution 
from  the  back  side  is  assumed  negligible,  so  that  only  a  consideration 
of  the  forces  acting  on  the  front  side  is  of  interest.  For  simplicity, 
the  pressure  forces  are  assumed  to  be  given  by  Newtonian  theory.  The 
Newtonian  force  coefficients  as  obtained  from  Reference  41  are: 
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where 


sin  6  =  tan  a/tan  a 
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The  center  of  pressure  for  Newtonian  theory  is  given  by  the  expression 
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These  expressions  are  also  applicable 
if  6  is  set  equal  to  n/2  for  a  <  a. 
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=  slant  length) 

to  the  small  angle  of  attack  case, 
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AC  large  angle  of  attack  It  will  be  assumed  that  Infinite  swept 
cylinder  theory^*-44  can  be  applied  to  the  slender  cones  of  Interest. 
Solutions  are  obtained  by  assuming  local  similarity  and  satisfying  the 
requirement  that  the  chordwise  and  spanwise  boundary  layer  development 
are  not  independent  for  variable  fluid  properties.  In  the  nomenclature 
of  Reference  44,  the  spanwise  and  chordwiae  skin  friction  shear  forces 
are,  respectively; 
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The  development  of  the  spanwise  skin  friction  solution  into  the 
axial  skin  friction  force  coefficient  is  as  follows:  The  definition  of 
the  pressure  gradient  parameter,  0,  is 
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For  an  isothermal  wall,  substitution  of  7.22  into  7.29  gives 
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The  axial  skin  friction  force  coefficient  is  given  by 
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This  integral  cannot  be  integrated  in  closed  form  unless  suitable  simplify¬ 
ing  distributions  for  the  local  velocity  and  density  are  assuned.  For  this 
reason  an  average  or  barred  value  of  p  will  be  assumed  and  a  linear  vari¬ 
ation  of  u„  with  ♦  such  that  Ug/U^,  =  2/n  *  sin  (ct  +  9o) .  Then 
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The  axial  skin  friction  force  coefficient  then  becomes 
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Solutions  for  g'w  as  a  function  of  P  and  heat  transfer  have  been  tabulated 
in  Reference  44.  To  the  accuracy  employed  in  the  previous  assumptions  for 
py,  and  ue ,  it  will  also  be  appropriate  to  assume  an  average  value  of  g'w 
such  that  a  linear  variation  of  p  with  ♦  gives 
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This  assumption  in  connection  with  Equation  6.3  gives  for  the  final  inte¬ 
gration: 
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The  only  undetermined  value  in  Equation  7.38  is  that  of  the 
average  density,  pw.  A  modified  Newtonian  pressure  distribution  with  0, 
and  a  linear  viscosity-temperature  relationship  gives 
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where  the  total  pressure  ^t0  is  the  value  obtained  at  the  stagnation 
line  from  the  normal  Mach  number  component.  Substitution  of  7.38  into 
7.37  and  algebraic  rearrangement  gives  for  the  final  solution: 


(1  -  tg)  sin  (a  +  Sg)  cos  (6g 


Re  sin  a 
"L 


7.40 

The  values  obtained  for  g'w  are  presented  in  Figure  7-1  for  a  Prandtl 
number  of  1.0. 


The  total  axial  force  coefficient  is  given  by 
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The  skin  friction  contribution  to  the  normal  force  coefficient 
may  be  obtained  in  a  similar  manner  to  that  for  the  axial  coefficient  by 
using  Equation  7.30  in  place  of  7.29. 
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Values  of  f"  for  Pr  =  1  are  given  in  Figure  7-2.  It  is  noted  that  the 
variation  ofWf"w  with  Prandtl  number  is  much  larger  than  in  the  case  of 
g'w,  with f "w  decreasing  as  the  Prandtl  number  decreases. 

The  moment  coefficient,  based  on  the  cone  axial  length  as  the 
reference  length  follows  from  Equation  7.42  by  multiplying  the  integrand 
by  the  moment  arm,  and  is  given  by 


CF  _  zl  3i_  7.44 
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The  center  of  pressure  is  again  given  by  Equation  7.24  where 
the  total  normal  and  moment  coefficients  are  obtained  by  summing  up  the 
individual  contributions.  The  moment  about  the  c.g.  location  is  given 
by 
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7.6  COMPARISON  OF  EXPERIMENT  AND  THEORY 

The  analytical  solutions  presented  have  been  compared  with 
experiment  in  Figures  7-3  to  7-10.  Also  shown  for  comparison  are  the 
Newtonian  estimates  of  the  aerodynamic  forces.  The  region  of  inter¬ 
mediate  angle  of  attack  are  estimated  by  a  smooth  curve  arbitrarily 
faired  between  the  small  and  large  angle  of  attack  solutions. 

Data  from  Reference  45  for  an  insulated  9-degree  half-angle 
cone  at  ~  5.8  are  used  for  the  experimental  data  in  Figures  7-3  to 
7-5.  These  data  are  approximated  by  Newtonian  theory  for  the  normal 


7-14 


7-15 


S12018 

COEFFICIENT) 


COEFFICIENT,  C  (CG)  AXIAL  FORCE  COEFFICIENT 


AXIAL  FORCE  COEFFICIENT 


FIGURE  7-6.  COMPARISON  OF  THEORY  AND  EXPERIMENT(AXIAL  FORCE  COEFFICIENT) 
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FIGURE  7-7.  COMPARISON  OF  THEORY  AND  EXPERIMENT (NORMAL  FORCE  COEFFICIENT) 
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and  moment  forces,  but  are  seriously  underestimated  by  Newtonian  theory  In 
the  case  of  the  axial  forces.  The  latter  condition  Is  caused  by  large 
viscous  Interaction  effects  at  small  a,  and  non-negliglble  skin  friction 
contributions  at  large  a.  According  to  Reference  45,  trailing  edge  drag 
effects  could  also  be  present  at  large  a  which  might  partially  account  for 
the  Increase  of  experimental  values  over  the  predicted  ones. 

The  contribution  of  the  back  side  of  the  cone  to  the  forces  was 
found  to  be  very  important  in  the  region  a  >  a  at  small  angle  of  attack. 

The  method  described  In  Section  7.3  was  employed  to  estimate  the  various 
contributions  to  the  total  force.  It  was  found  that  the  magnitude  of  the 
skin  friction  forces  on  the  back  side  was  of  the  same  order  as  that  on  the 
front  side.  The  viscous  interaction  effects  tend  to  compensate  for  the 
decrease  in  dynamic  pressure.  In  the  case  of  the  axial  force  (Figure7-4), 
this  provides  a  continuous  increase  in  CA  with  increasing  a,  corresponding 
to  the  experimental  results.  In  the  case  of  the  normal  and  moment  forces 
(Figures  7-3  and  7-5)  ,  the  skin  friction  forces  on  the  front  and  back  side 
tend  to  cancel  each  other,  thereby  contributing  very  little  to  the  net 
forces  in  the  low-density  flow.  An  examination  of  the  pressure  forces  acting 
on  the  cone  was  sufficient  for  an  approximate  analysis  of  the  normal  and  moment 
forces. 

The  unavailability  of  high  Mach  number  data  characteristic  of 
re-entry  has  prompted  Aeronutronic  to  initiate  a  series  of  tests  to  correct 
this  deficiency.  The  first  information  obtained  on  sharp  cones  of  6-1/3- 
and  9-degree  semi-vertex  angle  was  obtained  in  the  AEDC  Von  Karman  Gas 
Dynamics  Facility  in  Hotshot  Tunnel  Number  1  by  ARO ,  Inc.  Aerodynamic 
forces  were  measured  for  0  <  a  <  15  degrees.  Test  conditions  and  geometry 
were  as  follows: 
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This  data  has  been  examined  and  compared  with  theory  in  Figures  6-10  to 
6-14. 


In  Figure  7-6,  the  viscous  interaction  effects  become  apparent 
on  the  axial  force  coefficient.  Theory  predicts  a  crossover  of  the  6-1/3- 
and  9-degree  cone  curves,  with  the  6-1/3-degree  cone  having  greater  values 
of  CA  at  the  lower  Reynolds  numbers.  This  is  caused  by  the  increased 
skin  friction  due  to  the  transverse  curvature  effect.  Theory  also  predicts 
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a  change  In  the  crossover  to  higher  Reynolds  numbers  with  increasing  (2. 

The  data  were  obtained  in  this  region  where  the  crossover  occurs,  and  tend 
to  verify  the  predicted  trends.  In  the  case  of  the  normal  and  moment 
force  coefficients  (Figures  7-7  and  7-9) ,  the  data  for  the  9-degree  cone 
are  lower  than  the  Newtonian  estimation,  in  contrast  to  the  previous  data 
of  Reference  18,  but  in  agreement  with  prediction  at  small  at.  For  the 
6-1/3-degree  cone,  the  transverse  curvature  effect  is  sufficiently  strong 
to  increase  the  predicted  values  for  Re^.  =  10^  over  the  Newtonian.  The 
experimental  data,  however,  are  still  below  the  theoretical  estimates. 

As  the  angle  of  attack  continues  to  increase,  the  data  should  rise  above 
the  Newtonian,  similar  to  the  predicted  trend  in  Figure  7-3.  At  the  large 
angle  of  attack  where  the  Newtonian  estimate  adequately  accounts  for  the 
pressure  contribution  to  the  forces,  any  increase  in  the  data  over  the 
Newtonian  must  be  attributed  to  skin  friction  shear  forces.  These  should 
always  be  relatively  small  in  comparison  to  the  pressure  forces,  although 
not  necessarily  negligible. 

The  methods  outlined  previously  have  also  been  applied  to  con¬ 
ditions  obtained  in  the  Cornell  Aeronautical  Labs'  shock  tunnel  for  hyper¬ 
sonic  flow  data.  The  teste  were  conducted  for  Aeronutronic  under  BSD 
contract  to  ARO ,  Inc.,  as  an  extension  to  a  current  ARO  contract.  The 
comparison  of  the  theory  with  experiment  is  presented  in  Figure  7-11  for 
the  axial  force  coefficient,  and  in. Figures  7-12  ard  7-13  for  the  normal 
and  moment  force  data,  respectively.  Tabulated  data  are  given  in  Table  III 
for  all  of  the  force  measurements.  The  tabulated  data  also  include  measure¬ 
ments  on  a  blunt  6.34°  half-angle  cone. 

The  agreement  between  theory  and  experiment  for  the  sharp  cones 
is  found  to  be  good.  The  viscous  interaction  and  transverse  curvature 
effects  are  strongly  exhibited  in  the  axial  force  data  (Figure  7-11)  .  At 
the  same  length  Reynolds  number,  the  data  for  6.34°  cone  are  seen  to  cross¬ 
over  the  data  for  the  9-degree  cone  as  <2  increases,  and  actually  becomes 
larger  than  the  9-degree  cone  data.  The  zero  angle  of  attack  data  exhibit 
the  same  characteristics.  It  is  noted  that  the  disagreement  between 
Newtonian  theory  and  experiment  in  the  case  of  the  moment  force  data  is 
not  as  bad  as  might  seem  from  a  cursory  examination.  In  the  form  pre¬ 
sented,  the  accuracy  of  the  data  is  probably  ±  50  percent  or  greater,  and 
is  really  a  measure  of  the  agreement  between  theory  and  experiment  of  the 
center  of  pressure  location.  On  this  basis,  the  agreement  is  considered 
good . 


A  comparison  of  the  sharp  and  blunt  6.34°  cone  axial  force  data 
are  given  in  Figure  7-14.  The  viscous  interaction  effects  become  even 
more  evident  for  this  case,  as  the  axial  force  data  for  the  sharp  cone 
has  Increased  to  equal  that  of  the.  blunt  cone  over  the  Reynolds  number 
range  examined. 
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FIGURE  7-11.  AXIAL  FORCE  COEFFICENT  VERSUS  ANGLE  OF  ATTACK 
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FIGURE  7-12.  NORMAL  FORCE  COEFFICIENT  VERSUS  ANGLE  OF  ATTACK 
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TABLE  III 

l  ,  ■  Db 
ref  » 

A  -  ■  A_ 
ref  a 

i  ■  model  length  from  theoretical 


nose 


RUN  NO. 

0, 

DEGREES 

rn/rb 

Moo 

Rfioo 

FT'1  . 
x  10‘6 

a,  C 

DEGREES 

C  Cj,  BASE 

(about  DIAMETER 
2/3  i)  INCHES 

1 

9 

0.03 

14.86 

0.193 

0 

0.1228 

0.0008  -0.0002 

3.0 

2 

9 

0.03 

14.72 

0.180 

5 

0.1301 

0.1617  -0.0127 

3.0 

3 

9 

0.03 

14.70 

0.182 

10 

0 . 1485 

0.3165  -0.0291 

3.0 

4 

9 

0.03 

14.67 

0.172 

15 

0.1891 

0.5006  -0.0398 

3.0 

6 

9 

0  .03 

14.65 

0.171 

15 

0.1943 

0.4924  -0.0403 

3.0 

5 

9 

0.03 

14.67 

0.169 

25 

0.2756 

0.8821  -0.0555 

3.0 

7 

9 

0.03 

13.63 

0.022 

0 

0.2196 

0.0014  -0.0003 

3.0 

12 

9 

0.03 

13.64 

0.023 

p 

0.2363 

0.0087  0.0015 

3.0 

8 

9 

0.03 

13.70 

0.023 

10 

0.2604 

0.3077  -0.0363 

3.0 

9 

9 

0.03 

13.75 

0.024 

20 

0.3416 

0.6686  -0.0652 

3.0 

10 

9 

0.03 

13.73 

0.024 

30 

0.4888 

1.2528  -0.0966 

3.0 

11 

9 

0.03 

13.87 

0.024 

40 

0.5758 

1.7534  -0.1509 

3.0 

13 

9 

0.03 

14.36 

0.068 

0 

0.1482 

-0.003  0.0002 

3.0 

14 

9 

0.03 

14.30 

0.065 

15 

0.2373 

0.5151  -0.0523 

3.0 

15 

9 

0.03 

14.31 

0.065 

40 

0.4456 

1.6165  -0.1263 

3.0 

16 

6.34 

0.03 

14.66 

0.172 

0 

0.0988 

-0.0017  0.009 

2.375 

21 

6.34 

0.03 

14.64 

0.168 

0 

0.1031 

-0.0008  0.0005 

2.375 

17 

6.34 

0.03 

14.65 

0.171 

5 

0.1177 

0.1672  -0.0066 

2.375 

It 

£  .34 

0.03 

14.64 

0.170 

10  • 

0.1457 

0.3480  -0.0117 

2.375 

19 

6.34 

0.03 

14.61 

0.167 

15 

0.1917 

0.5960  -0.0269 

2.375 

22 

6.34 

0.03 

14.63 

0.168 

15 

0.1886 

0.5638  -0.0205 

2.375 
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TABLE  III  (Continued) 


ref  "B 
Aref  “  *B 

i  -  model  length  from 
theoretical  nose 


RUN  NO. 

o, 

DEGREES 

VS 

Moo 

Reoo 
FT'1 
x  10'6 

a*  CA 

DEGREES 

C»  .&>„« 
2/3i) 

BASE 

DIAMETER 

INCHES 

'*20 

| 

6.34 

0.03 

14.64 

0.171 

25 

0.2828 

1.0631  -0.0232 

2.375 

23 

6.34 

0.03 

13.79 

0  .024 

0 

0.2323 

0.0106  0.0087 

2.375 

27 

6.34 

0.03 

14.36 

0.067 

0 

0.1413 

-0.0003  -0.0046 

2.375 

28 

6.34 

0.03 

13.65 

0.022 

0 

0.2219 

0.0065  0.0020 

2.375 

24 

6.34 

0.03 

13.72 

0.023 

10 

0.2720 

0.3460  -0.0294 

2.375 

25 

6.34 

0.03 

13.65 

0.022 

20 

0.4332 

0.9477  -0.0489 

2.375 

30 

6.34 

0.03 

13.77 

0.024 

20 

0.4239 

0.8330  -0.0297 

2.375 

*26 

6.34 

0.03 

13.69 

0.023 

30 

0.5516 

1.5918  -0.0443 

2.375 

*29 

6.34 

0.03 

13.72 

0.023 

40 

0.6720 

2.3421  -0.0381 

2.375 

*31 

6.34 

0.03 

14.59 

0.171 

35 

0.3756 

1.6705  -0.0177 

2.375 

*32 

6.34 

0  .03 

14.54 

0.162 

45 

0.4529 

2.3402  -0.0053 

2.375 

33 

6.34 

0.144 

14.59 

0.162 

0 

0.1007 

-0.0030  -0.0007 

2.375 

34 

6.34 

0.144 

14.59 

0.166 

5 

0.1240 

0.1549  -0.0384 

2.375 

45 

6.34 

0.144 

14.62 

0.165 

5 

0.1275 

0.1548  -0.0389 

2.375 

35 

6.34 

0.144 

14.60 

0.170 

10 

0.1531 

0.3256  -0 .0426 

2.375 

36 

6.34 

0.144 

14.61 

0.166 

20 

0.2432 

0.7752  -0.0522 

2.375 

46 

6.34 

0.144 

14.63 

0.166 

40 

0.4115 

2.1096  -0.0546 

2.375 

37 

6.34 

0.144 

13.63 

0  .023 

0 

0.2217 

0.0015  0.0040 

2.375 

38 

6.34 

0.144 

13.72 

0.024 

5 

0.2447 

0.1689  -0.0426 

2.375 

39 

6.34 

0.144 

13.67 

0.024 

10 

0.2928 

0.3694  -0.0535 

2.375 

40 

6.34 

0.144 

13.67 

0.023 

20 

0.4194 

0.8791  -0.0846 

2.375 

42 

6.34 

0.144 

13.80 

0.024 

30 

0.5567 

1.6587  -0.1105 

2.375 

41 

6.34 

0.144 

13.76 

0  .024 

40 

0.6576 

2.4687  -0.1669 

2.375 

*  Validity  of  moment  force  data  is  questionable. 
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FIGURE  7-14 .  COMPARISON  OF  C *  FOR  SHARP  AND  BLUNT  SLENDER  CONES 
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In  conclusion,  a  reasonable  engineering  approximation  to  the 
forces  acting  on  sharp,  slender  cones,  even  with  low-density  flow  effects, 
would  be  to  apply  Newtonian  estimates  for  the  normal  and  moment  forces  in 
hypersonic  flow.  The  axial  forces,  however,  must  consider  the  viscous 
interaction  effects  which  greatly  increase  both  the  pressure  and  skin 
friction  contributions. 
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CHAPTER  8 
NOMENCLATURE 


al*  a2’  a3 

V 


A 

\ 

A,  B 
C 


C 


M 


d,  G 


D 

E 


coefficients  given  by  Equation  5.20 

parameters  describing  shock  shape  in  hypersonic  small 
disturbance  theory  (equation  3.26) 

area 

cone  base  area 

coefficients  given  in  Equation  2.2 

coefficient  in  linear  viscosity-temperature  relationship, 
M  /  Me  =  C  T/Tg 

specific  heat 

pressure  coefficient 

local  skin  friction  coefficients  shear  force/qmA 

c 

mean  skin  friction  coefficient  =  total  shear  force/q  A 

00  Q 

drag  coefficient,  P/q^  A^ 
axial  force  coefficient 
normal  force  coefficient 
moment  force  coefficient 

coefficients  expressing  heat  transfer  effects  in  viscous 
interaction  theories,  Equations  2.2,  2.7 

cone  base  diameter 

Equation  5  .9 
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NOMENCLATURE  (Continued) 
coefficient  given  by  Equation  2.22 
Equations  3.13,  3.14 

function  describing  Induced  pressures  in  strong 
interaction  region  (equation  2.8) 

function  describing  boundary  layer  thickness  in  strong 
interaction  region  (equations  2.2a,  2.2b) 


h 


static  enthalpy 

also,  heat  transfer  coefficient 


L 

f 


total  enthalpy 

integral  given  in  Equation  2.9 

coefficients  given  in  Equations  2.3  and  2.4 

hypersonic  similarity  parameter,  Mo,^ 

also,  constant  of  proportionality  in  boundary  layer 

growth,  6*  *  Kx11 

slant  length 

length 


M 


Mach  number 


m 

n 


exponent  describing  power  law  body,  r^  a  x1" 
exponent  describing  boundary  layer  growth,  5*  dx" 


P 


pressure 


Pr 


Prandtl  number 


Q/A 


heat  flux 


q 


dynamic  pressure 
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R 

Re 

r 

St 

T 

T* 

u,  v 
x,  y 

z 

a 

0 

7 

6 

6* 

e 

0 

0 

s 

9b 

X 


NOMENCLATURE  (Continued) 

■  cylindrical  radius  of  shock  wave 
also,  gas  constant 

also,  equivalent  radius  of  curvature 
=  Reynolds  number 

■  local  cylindrical  radius 

■  Stanton  number 


■  temperature 
*  reference  temperature 


*  velocity  components 

*  coordinate  system  dimensions,  x  taken  along  slant  length, 
and  y  normal  to  surface 

«  compressibility  factor  in  equation  of  state 


*  angle  of  attack 


35  shock  shape  parameter  —  -  1 

i  ,  ,  .  m 

also,  pressure  gradient  parameter, 


Izl  • 

7 


n 

1+n 


*  ratio  of  specific  heats 

■  boundary  layer  thickness 

■  displacement  thickness 
»  density  ratio,  P\/p  2 

■  total  flow  deflection  angle 

■  angle  between  local  tangent  to  shock  wave  and  body  axis 
(shock  angle) 

«  local  body  angle 

■  non-dimensional  distance,  Equation  3.6 
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NOMENCLATURE  (Continued) 

M 

- 

viscosity 

5 

SB 

Equation  3.15 

P 

a 

density 

a 

as 

semi-vertex  angle  of  sharp  cone 

X 

m 

viscous  interaction  parameter,  M  </"c7  vRe 

0 

m 

circumferential  body  angle  (0"  0  at  stagnation  line) 

A 

m 

angle  of  yaw 

A 

m 

shock  standoff  distance 

Subscripts 

0,  orig 

SB 

original  inviscid  value,  also 
stagnation  value 

w 

a 

wall 

c 

a 

cone 

t 

a 

total,  or  stagnation  value 

e 

a 

local  value  at  edge  of  boundary  layer 

1,2 

a 

values  before  and  after  shock  wave,  respectively 

00 

a 

freestream 

b 

a 

body,  or  base 

M 

a 

Mangier  solution 

Aw 

a 

adiabatic  wall 

P 

- 

local  value 

L 

a 

slant  length 
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NOMENCLATURE  (Continued) 

Subscripts  (continued) 


M 

- 

normal 

s 

■1 

shock 

ref 

- 

reference  value 

FM 

m 

free  molecule  flow 

P 

- 

pressure  contribution 

F 

. 

skin  friction  contribution 
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APPENDIX 

PROCEDURE  FOR  CALCULATING  SHOCK  SHAPE  AT  LARGE  ANGLES  OF  ATTACK 


The  slender  cone  Is  assumed  to  act  as  a  cylinder  with  an  equiva¬ 
lent  radius  of  curvature,  Rb  (equation  6.3  of  text).  At  large  angles  of 
attack  a  swept  cylinder  is  approximately  equivalent  to  another  cylinder 
which  is  unswept  to  the  normal  Mach  number  component,  Mjj.  If  a  slender 
cone  is  considered  whose  equivalent  freestream  Mach  number  is  Mr 
■  Mx  sin  (@8  +  Cl)  ,  the  results  of  Love  may  be  applied  to  determine  the 
local  shock  slope  relative  to  the  direction  M^: 


**N 


also, 


Y/D 


SHOCK  WAVE 


SONIC  POINT  ON  SHOCK 


FIG.  A-l 


FIG.  A-2 


F/D 


A-l 
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The  input  quantities  required  to  compute  the  shock  wave  are: 
Mj,  -  Mm  sin  (0.  +a) 


€  ■  shock  angle  (with  respect  to  x  *  axis)  that  gives 

sonic  velocity  behind  the  shock  wave  (to  be  obtained 
from  NACA  Report  1135) 

Sdet  *  semi-vertex  angle  of  wedge  that  just  produces  shock 
detachm.nl  (NACA  Report  Number  1135) 

Cc  **  values  from  NACA  TN  4170  for  2-dimensional  circular 
1  J  face  (e.g.  cylinder  normal  to  flow),  see  following 

tables: 


“n 

C 

c 

q  ,  degrees 

1.0 

... 

0 

1.5 

— 

4 

1.6 

0.780 

— 

2.0 

.806 

11.6 

2.4 

.816 

— 

2.5 

— 

20.0 

2.8 

.823 

— 

3.0 

— 

27.5 

3.5 

— 

34.5 

3.6 

.836 

— 

4.0 

— 

41.0 

4.4 

.845 

— 

4.5 

— 

46.5 

5.0 

— 

52.0 

5.2 

.850 

— 

5.5 

— 

57.0 

6.0 

.854 

62.0 

7.2 

.856 

— 

OO 

0.857 

-137. 
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The  equations  for  the  shock  shape  are  as  follows: 


77  -  0.50  C  cot  6.  _ 

a  c  det 


o 

d' 


P  yp2  tan2g8  -  1  ( JT  + 

2  /  2  2  1 
°>  tan  e  -P  w  P  tan  c  -  1  +  tan  n 

8*8 


Y_ 

d' 


1 

P 


m2  -W 


D  *  d7  C0S&det 


F 

D 


>  -  d^|  COS&det 


sin2&  +  C  cos2  6 

det  c  det 


2  sin  6 


det 


The  local  shock  slope  is 


A.  1 


A. 2 


A. 3 


A. 4 


A. 5 


±L  = 

dx 


A. 6 


Experimental  data  indicates  that  the  shock  detachment  distance  given  by 
equation  6.7  should  be  used.  Then 


d’/2  -  R,  6  , 

b  cos  det 


and  in  place  of  equation  A.l,  we  have 


*1 

d' 


+  (1  -  sin  5 


A. 7 


A. la 
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